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Section  1 


INTRODUCTION 

This  report  describes  the  results  of  theoretical  and  experimental  investigations 
of  the  coherent  anti-Stokes  Raman  scattering  (CARS)  process  of  molecular  gases 
which  were  conducted  under  USAF  Contract  F33615-77-C-2019. 

A  program  is  in  progress  at  the  AFWAL  Aero  Propulsion  Laboratory  (APL)  to 
evaluate  combustor  models  in  environments  which  simulate  many  of  the  features 
in  gas-turbine  combustors.  This  program  involves  the  selection  and  evaluation 
of  appropriate  diagnostic  techniques  for  making  time-averaged  and  time-resolved 
point  measurements  of  velocity,  temperature,  and  major  species  concentration 
in  simulated,  practical  combustion  environments. 

Of  the  various  optical  diagnostic  techniques,  coherent  anti-Stokes  Raman 
spectroscopy  (CARS)  is  considered  to  be  one  of  the  more  promising  methods  for 
making  temperature  and  species-concentration  measurements  in  combustion  envi¬ 
ronments.  Favorable  results  have  been  obtained  by  various  investigators  using 

1-4 

laboratory-type  flame  sources.  In  addition,  a  preliminary  evaluation  of 

this  technique  using  the  large-scale  AFWAL/APL  combustor  has  been  made,^’^  and 

the  results  are  very  encouraging  in  view  of  the  large  scale,  hostile,  gas- 

turbine-type  combustion  environment  in  which  the  sensitive  optical  system  was 

operated.  Researchers  at  United  Technologies  Research  Center  have  made 

successful  temperature  measurements  using  the  BOXCARS  technique  in  a  large- 

scale  furnace-type  combustor.^  Investigators  at  Shell  Research  Limited 

Thornton  Research  Centre  have  demonstrated  the  feasibility  of  using  CARS  to 

8 

measure  temperatures  in  a  reciprocating  engine. 

The  overall  objective  of  this  contract  was  to  establish  the  CARS  technique  as 
a  valid  combustion  diagnostic  tool  through  fundamental  laboratory  studies  and 
to  demonstrate  the  practicality  of  CARS  through  in  situ  testing  of  an  environ¬ 
mentally  hardened  CARS  system  developed  during  the  program. 
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Section  2  contains  a  short  comprehensive  review  of  the  CARS  theory  to  serve 
primarily  as  background  material  for  the  remaining  sections. 

Section  3  describes  the  fundamental  laboratory  equipment  and  studies  which 
were  undertaken.  These  include:  1)  linewidth  determinations  by  integrated 
power  measurements;  2)  flame-breakdown  measurements;  3)  single-shot  thermometry; 
A)  reference-scheme  studies,  5)  simultaneous  temperature  and  species- 
concentration  measurements,  and  6)  turbulence  effects  and  an  experimental 
arrangement  for  their  compensation. 

Section  A  describes  the  construction,  evaluation,  and  applications  of  the 
hardened  CARS  system  for  combustion  diagnostics  on  the  AFWAL/APL  large-scale 
practical  combustor.  This  includes  temperature  and  species-concentration 
profiles  of  the  combustor,  comparisons  of  the  CARS  results  with  probe  measure¬ 
ments,  and  an  overall  evaluation  of  the  performance  of  the  system  in  a  hostile, 
dynamic  environment. 
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Section  2 


THEORY 

In  this  section,  highlights  in  the  theory  of  the  CARS  process  are  presented  in 
order  to  bring  out  the  essential  features  in  the  analyses  of  the  experimental 
data.  A  more  complete  description  of  the  theory  can  be  found  in  Refs.  9-10, 

Due  to  the  nonlinear  nature  of  the  response  of  a  medium  to  strong  external 
stimuli,  the  polarization  induced  in  a  medium  subjected  to  intense  light  waves 
contains  not  only  the  frequency  components  of  the  incident  waves  but  also 
linear  combinations  of  these  frequencies.  This  induced  polarization  containing 
many  frequency  components  acts,  in  turn,  as  the  source  of  a  secondary  wave 
containing  the  same  number  of  frequency  components. 


In  CARS  experiments,  two  laser  beams  are  used  as  the  incident  waves,  one 
wave  having  the  higher  frequency  at  (called  the  pump  beam)  and  the  other 
having  the  lower  frequency  (called  the  Stokes  beam).  Of  the  many  secondary 
waves  generated  in  the  medium,  the  one  used  as  the  signal  in  a  CARS  experiment 
has  the  frequency  at  2 -  ui^  (the  anti-Stokes  beam).  As  a  result  of  a 
resonance  which  occurs  in  the  induced  polarization  when  the  difference 
frequency  of  the  two  laser  beams  is  equal  to  a  Raman-active  vibrational 
frequency  (u^)  of  a  molecule  in  the  medium  (  i.e.,  u>^  -  the  CARS 

signal  generated  at  =  2u>^  -  (=  +  uO  undergoes  a  resonant  enhancement 

in  intensity.  Depending  upon  the  linewidth  of  the  laser  beams  employed  in 
the  experiment,  slightly  different  versions  of  the  theory  are  applicable  in 
the  analysis  of  the  CARS  process. 


In  the  case  of  monochromatic  incident  waves,  the  nonlinear  polarization 
induced  at  the  anti-Stokes  frequency  is 


P<3)(t)  =  3X(3)eJ  E;  e^(2kl  "  V  -  C^-^tl 


(1) 


with 


(3) 
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1 _ 

(“l  ~  w2^ 
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and 


NA  T  =  N  -  N  , 
v,J  v,J  v+l,J 


(3) 


where  is  the  third-order  CARS  susceptibility;  is  the  nonresonant 

third-order  susceptibility;  do/dft  is  the  spontaneous  Raman  cross  section; 

ui  is  the  angular  frequency  of  the  Q-branch  Raman  transition;  N  is  the 
v ,  J 

molecular  number  density;  and  subscripts  v  and  J  indicate,  respectively, 
vibrational  and  rotational  quantum  numbers  of  the  energy  level  to  which 
various  parameters  are  related.  Moreover,  in  thermal  equilibrium,  the 
molecular  distribution  is  Boltzmann,  i.e.. 


N 


v,  J 


N  (2J  +  1)  gT  -Ev  j/kT 
jj  e 


(4) 


where  Q  is  the  partition  function  and  is  the  nuclear-spin  degeneracy. 


When  Maxwell's  equations  which  describe  the  propagation  of  electromagnetic 
waves  are  solved  for  the  anti-Stokes  wave  after  substitution  of  the  third- 
order  polarization  as  given  in  Eq.  (1),  it  can  be  shown  that  the  intensity 


of  the  anti-Stokes  wave  generated  in  an  interaction  length  x  is  given,  in  the 
case  of  a  phase-matched  interaction,  by 


X3  =l 


,  2  x  2 
4t,  u>3\ 


|3X0)|2  I2  I/ 


(5) 
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Since  the  CARS  intensity  is  proportional  to  1^  I^,  one  can  achieve  a  degree 
of  spatial  resolution  in  the  measurement  by  focusing  the  incident  laser  beams. 
The  lateral  dimension  of  the  laser  beams  (called  spot  size  for  Gaussian  beams) 
gives  the  lateral  resolution,  and  the  length  of  the  focal  zone  (confocal 
parameter  for  Gaussian  beams)  gives  the  longitudinal  resolution  along  the 
beam  direction.  In  such  situations  it  is  more  appropriate  and  convenient  to 
deal  with  the  total  power  of  the  anti-Stokes  beam  rather  than  the  intensity. 
Thus,  in  the  case  of  the  focused  Gaussian  beams,  the  anti-Stokes  power  generated 
can  be  approximated  by 
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Since  the  measurement  of  the  power  of  the  anti-Stokes  wave  as  given  by 

(3) 

Eq.  (6)  will  yield  information  for  x  which,  in  turn,  is  related  to 
N  by  Eqs.  (2)  -  (4),  it  should  be  possible  to  determine  the  number  density 
for  various  states  of  any  Raman-active  molecular  species  present  in  the 
combustion  media.  This,  in  fact,  forms  the  basis  for  the  application  of  the 
CARS  process  to  gas  diagnostics — particularly  in  combustion  media. 


In  the  case  of  nonmonochromatic  incident  waves,  the  analysis  is  best  carried 
out  in  the  frequency  domain.  Since  in  actual  experimental  systems  the 
incident  beams  are  not  really  monochromatic,  this  case  is  more  representative 
of  a  real  experimental  situation.  In  fact,  the  data  presented  in  this  paper 
were  obtained  using  Gaussian  beams  from  a  single-mode  ruby  laser  at  the 
frequency  oi  and  a  broadband  dye  laser  centered  at  frequency  u>  .  The  measured 
linewidth  of  the  single-mode  pump  laser  was  -  0.034  cm  ;  consequently,  it 
was  assumed  to  be  monochromatic.  This  approximation  turns  out  to  be  satisfac¬ 
tory  for  the  applications  under  discussion.  In  this  approximation,  the 
spectral  power  density  of  the  anti-Stokes  wave  becomes 


(uj^  +  6)  =1 


4uu)^u)^ 


3X(3V 


+  6)|  2  P2  P2(~u>2  +  6) 


Equation  (7),  when  convolved  with  an  instrument  slit  function,  then  can  be 
used  to  calculate  normalized  CARS  spectra  as  a  function  of  temperature. 


For  number-density  measurements  the  bandpass  of  the  signal-analyzing  spec¬ 
trometer  was  increased  to  allow  measurement  of  the  energy  in  the  entire 
Q-branch.  If  the  spectral  power  density  of  the  broadband  dye  laser  is 
assumed  to  be  constant  over  the  Q-branch,  then  the  integrated  anti-Stokes 
power  (PINT)  becomes 
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(8) 
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In  evaluating  the  integral  in  Eq.  (8),  the  Raman  linewidth  (F)  was  assumed 
to  be  constant  for  all  lines  and  the  contribution  of  the  hot  bands  was 
neglected. 


Species  number  densities  were  determined  by  comparing  P^T  measured  in  the 
combustion  gas  with  (PT„„)  values  generated  in  a  calibration  cell  containing 
a  known  number  density  of  the  same  species.  The  calculation  of  F(T)  required 
knowledge  of  the  combustion  gas  temperature,  and  F(Tc)  was  determined  from 
the  temperature  of  the  calibration  gas  (Tc).  In  addition,  knowledge  of 
as  a  function  of  pressure,  temperature, and  composition  was  required.  Since 
those  data  were  not  available,  only  an  estimated  value  of  T  could  be  used. 
This  value,  obtained  by  assuming  the  linewidth  to  be  proportional  to  the 
collisional  frequency  of  the  molecule,  is 

r  «  a  v(T)  N  n(T)  (9) 

eff  o  i 

where  a  is  assumed  to  be  a  constant  optical  cross  section;  v(T)  is  the  average 

o 

molecular  velocity;  is  the  total  gas  density;  and  n(T)  is  a  broadening 
parameter.  From  laboratory  measurements  the  value  of  n(300  K)  for  and  09 
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was  determined  to  be  ~  0.8.  At  combustion  temperatures  the  value  of  n  was 

assumed  to  be  equal  to  1  due  to  the  increased  separation  between  populated 

J  levels.  By  incorporating  these  assumptions  for  into  Eq.  (8)  and 

normalizing  by  taking  the  ratio  P,„_,/(PT„m)  ,  one  obtains  the  expression  used 

ini  INI  c 

for  determining  species  number  density 


N 


fc.  (hm—\  (im 

kTc  V'WcJ  \F(T) 


v3-  V* 

k-2T'5pJ 


(10) 


where  p  and  pc  are  the  pressures  of  the  combustion  and  calibration  gases, 
respectively,  and  k  is  the  Boltzmann  constant.  The  quantity  indicated  by  the 
second  bracket  in  Eq.  (10)  can  be  viewed  as  a  correction  factor  to  the 
number  density  which  is  required  to  offset  the  effects  upon  F(T)  and  T  .. 
between  Tc  and  T.  This  factor  is  plotted  in  Fig.  1  for  and  0^  at 
T  =  300  K  and  p  =  p  =1  atm.  The  number  densities  obtained  from  this 
equation  were  then  converted  to  mole  fractions  using  the  ideal  gas  law  and 
the  CARS-determined  flame  temperature. 


7 


NUMBE 


Figure  1.  Plot  of  Number-Density  Correction  Factors 
as  a  Function  of  Temperature. 
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Section  3 


FUNDAMENTAL  LABORATORY  WORK 


3 . 1  EXPERIMENTAL  APPARATUS 

Ruby-Based  CARS  System 

The  ruby-based  CARS  system  is  shown  schematically  in  Fig.  2.  The  system 
consists  of  a  ruby-laser  oscillator /amplifier ,  a  ruby-pumped  infrared  dye 
laser,  monochromator/detector  and  electronics,  and  various  optics  for  guiding 
the  beams.  In  addition,  a  He-Ne  laser  beam  is  used  for  aligning  the  entire 
CARS  system.  Basically,  the  system  operates  as  follows:  A  single-mode 
(transverse  and  longitudinal)  ruby-laser  beam  emanating  from  the  ruby 
oscillator  is  amplified  in  the  ruby  amplifier  and  split  into  two  parts  by 
means  of  a  beamsplitter  (BS1) .  One  part  (70%)  of  the  ruby  beam  is  weakly 
focused  and  directed  into  the  dye  cell  slightly  off-axis  for  pumping  the  dye 
laser.  The  dye-laser  output  is  combined  collinearly  with  the  other  part  of 
the  ruby-laser  beam  and  focused  into  a  reference  cell  containing  a  non¬ 
resonant  gas  under  high  pressure  (typically  30  atm).  This  collinear  super¬ 
position  of  the  pump  and  the  Stokes  beams  satisfies  the  phase-matching 
conditions  for  coherent  generation  of  anti-Stokes  radiation  in  gases.  The 
anti-Stokes  signal  generated  in  the  reference  cell  is  removed  from  the  laser 
beams  by  means  of  a  dichroic  mirror  (DM2),  and  the  laser  beams  are  refocused 
into  the  sample  cell  containing  the  medium  under  investigation.  The  ruby- 
laser  power  and  pulse  shape  are  continuously  monitored  by  means  of  a  PIN 
diode  detector  and  a  transient  digitizer  (Tektronix  R  7912 — not  shown  in  ‘rhe 
figure).  The  anti-Stokes  signals — one  from  the  reference  cell  and  the  other 
from  the  sample  cell — are  detected  and  processed  using  either  of  the  detection 
systems  shown  in  Fig.  3,  depending  upon  the  particular  experimental  require¬ 
ments.  These  two  detection  systems  are  discussed  below  along  with  their 
respective  modes  of  operation. 

The  broadband  dye  mode  of  operation  is  used  for  instantaneous  generation  of 
an  entire  Q-branch  spectrum  of  anti-Stokes  radiation  in  a  light  molecule 
such  as  and  using  a  single  laser  pulse.  The  dye  laser  is  (typically 


RUBY  USER 


RUBY  AMPLIFIER 


Figure  2.  Schematic  Diagram  of  the  CARS  System.  The  components  are  as 

follows:  Mirrors  Ml-48;  achromatic  lenses  L2-L4;  simple  lenses 
LI,  L5-L6;  filters  F1-F4;  dichroic  mirrors  DM1-DM2;  beam 
splitter  (7OR/30T)  BS1;  prism  PI;  Q-switch  Ql;  output  etalon  El; 
optical  multichannel  analyzer  OMA. 
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Signal-Processing  Electronics  for  CARS 


100  wave  numbers  wide),  and  the  generated  anti-Stokes  signals  from  the  sample 
and  reference  cells  are  dispersed  within  a  single  monochromator  and  detected 
by  a  dual-spectrum  optical  multichannel  analyzer  (OMA) .  Separation  of  the  two 
optical  signals  is  accomplished  by  focusing  the  anti-Stokes  beams  on  two 
different  points  on  the  entrance  slit  of  the  monochromator.  Because  of  the 
instantaneous  nature  of  the  generation  of  the  entire  spectrum,  this  mode  of 
operation  is  more  desirable  for  real-time  measurement  of  gas  temperatures. 

The  narrow-band  dye  mode  of  operation  is  employed  for  generation  of  a  point  in 
the  anti-Stokes  spectrum  of  a  gas  utilizing  a  narrow-band  dye-laser  line 
(typically  <  0.2  cm  The  anti-Stokes  signals  are  dispersed  in  two  separate 

monochromators  and  detected  by  photomultipliers.  The  photomultiplier-output 
pulses  are  integrated  in  an  RC  circuit  and  amplified  before  being  stored  in  a 
transient  recorder.  The  stored  pulses  subsequently  are  displayed  on  an 
oscilloscope  for  monitoring.  Since  the  spectral  power  density  of  the  narrow- 
band  dye  laser  generally  is  much  higher  than  that  of  the  broadband  laser,  the 
anti-Stokes  intensity  generated  is  higher  with  narrow-band  operation.  Further¬ 
more,  since  photomultipliers  are  used  in  place  of  the  OMA  as  shown  in  Fig.  2, 
the  overall  sensitivity  also  is  much  higher.  However,  in  order  to  generate 
the  spectrum  of  an  entire  band  in  this  mode,  hundreds  of  laser  pulses  scanning 
the  Stokes  region  are  required.  For  these  reasons,  this  mode  of  operation  is 
better  suited  to  species-concentration  measurements. 

Nd:YAG- Based  CARS  System 

Well  into  the  middle  of  this  program  of  study,  it  was  decided  that  the 
Nd:YAG-based  CARS  system  would  be  the  optimum  system  for  a  combustion  tool 
because  of  its  larger  power  levels  and  the  reliability  of  the  laser.  In  an 
effort  to  reduce  mechanical  vibrations  and  beam  movement,  the  laboratory 
Nd:YAG  CARS  system  was  incorporated  into  a  single  4  by  8  ft.  NRC  floating 
optical  table.  Figures  4(a)  -  4(c)  depict  the  resulting  arrangement.  The 
system  consists  of  a  Quanta-Ray  Nd:YAG  laser  whose  frequency-doubled  output  is 
used  to  pump  a  broadband  dye  laser  and  subsequently  is  used  in  the  CARS 
process.  The  broadband  dye  laser  consists  of  an  oscillator-amplifier 
combination  with  a  resulting  30%  lasing  efficiency.  An  optical  delay  line 
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figure  4(aJ.  Present  CARS  Experimental  Arrangement  Depicting  the  Sample 


Iris 


(M8,  M9,  M10)  is  incorporated  into  the  scheme  to  provide  a  temporal  match 
between  the  Nd:YAG-  and  dye-laser  pulses.  The  Nd:YAG  and  dye  beams  in  the 
BOXCARS  orientation  are  crossed  and  focused  by  lens  L5  into  the  premixed 
burner  flame.  Beam  splitter  BS5  splits  off  -  30%  of  the  three  beams  and 
forms  the  reference  leg.  This  path,  as  discussed  later,  is  used  for  normaliz¬ 
ation  and  species-concentration  determinations.  The  resulting  CARS  signals 
in  both  the  sample  and  reference  are  then  predispersed  by  the  Pellen  Broka 
prisms  and  coupled  into  the  SPEX  1702  monochromator  by  periscopes  PS2  and 
PS3 . 

The  second  major  modification  was  the  replacement  of  the  PAR  1203  OMA 
(vidicon)  with  a  Tracor-Northern  TN1710  DARSS  (Diode  Array  Rapid  Scanning 
Spectrometer)  shown  in  Fig.  5.  The  DARSS  unit  consists  of  a  linear  (one¬ 
dimensional)  diode  array  (reticon)  detector  element  in  conjunction  with  a 
photo-cathode  intensifier.  The  overall  sensitivity  of  the  device  is  similar 
to  the  OMA  unit,  which  is  about  50%  that  of  a  common  S-20  response  tube. 

This  was  found  to  be  more  than  adequate  for  these  combustion  studies.  The 
major  advantages  of  the  DARSS  over  the  OMA  are  the  larger  linear-dynamic 
range,  the  lack  of  blooming  or  cross-talk  between  channels  of  the  detector, 
and  the  better  electronic  support  offered  by  the  TN1710  control  unit. 

In  order  to  circumvent  the  single-dimensionality  of  the  Tracor-Northern 
DARSS  detector,  the  entrance  slit  of  the  monochromator  was  rotated  90°  from 
its  usual  vertical  position  into  a  horizontal  position.  This  allowed  the 
two  CARS  signals  to  be  aligned  on  different  regions  of  the  DARSS  detector. 

The  two  CARS  signals  were  then  recorded  simultaneously  by  the  DARSS  and 
analyzed  at  a  later  date.  This  configuration  allows  simultaneous  determin¬ 
ation  of  temperature  and  species  concentration  as  discussed  later. 

Preliminary  studies^  indicated  the  best  possible  experimental  scheme  for 
normalization  of  the  sample  and  reference  to  be  that  shown  in  Fig.  4. 

Temporal  studies  indicated  that  single-mode  operation  of  the  Nd:YAG  and 
dye  lasers  was  also  highly  desirable.  This  requires  the  use  of  Quanta-Ray's 
electronic-line-narrowing  device.  This  device  acts  to  force  the  electro- 
optical  Q-switch  of  the  Nd:YAG  to  behave  as  a  lossy  Q-switch,  which  permits 
a  slow  build-up  of  the  laser  pulse  over  a  relatively  long  period  of  time. 
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This  slow  build-up  allows  the  intracavity  elaton  to  be  more  efficient  in 
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discriminating  against  unwanted  inodes.  However,  the  free-running  nature 
of  the  cavity-pulse  build-up  of  this  device  presents  a  major  problem  in  the 
timing  of  the  detection  electronics. 

To  circumvent  the  possible  skew  and  lag  problems  which  this  caused,  the 
firing  of  the  lamps  of  the  Nd:YAG  was  synchronized  with  the  integration/ 
digitization  cycle  of  the  TN  1710-DARSS.  Figure  6  depicts  the  time  events 
which  occur  in  the  present  configuration.  The  DARSS  digitization/ integration 
timing  cycle  is  used  to  key  the  firing  of  the  Nd:YAG  flashlamps.  After  a 
period  of  50  -  300  ps,  the  Q-switch  opens  and  70  ns  later  the  laser  output 
occurs.  This  process  continues  at  the  hertz  rate  set  by  the  integration 
time  of  the  DARSS  cycle,  which  is  usually  set  for  10-Hz  operation.  When 
data  are  to  be  taken,  the  DARSS  is  reset  for  trigger  operation,  and  the  data 
control,  data  acquisition,  and  transient  digitizer  are  triggered.  The 
transient  digitizer  is  used  for  examining  the  mode  quality  of  the  Nd:YAG  to 
ensure  that  a  single  longitudinal  mode  will  be  obtained.  This  is  necessary 
since  the  ELN  device  produces  single-mode  operation  only  50%  of  the  time. 

If  the  pulse  is  indeed  single  mode,  the  data  are  accepted  and  the  process 
continues.  This  method  assures  that  all  data  are  taken  under  single-mode 
operation.  • 

Premixed  Burner  Arrangement 

The  burner  consists  of  a  premixed  Perkin-Elmer  assembly  No.  290-0107  and  an 
adjustable  nebulizer  No.  303-0352.  The  assembly  was  fitted  with  a  water- 
cooled  argon  sheath  which  was  employed  for  stabilizing  the  flame.  Both 
acetylene-air  and  propane-air  premixes  were  used  in  the  burner  for  flame 
studies.  The  operating  conditions  of  the  burner  for  both  fuels  are  listed 
in  Table  1.  The  assembly  was  mounted  on  two  translation  stages  which  allowed 
height  and  axial  profiles  of  the  flame.  The  observed  flame  consisted  of  an 
inner  and  an  outer  portion,  both  of  which  could  be  independently  controlled. 
Figure  7  depicts  the  flame  that  was  obtained  with  the  burner  when  acetylene 
was  used  as  the  fuel.  The  overall  height  of  the  flame  was  greatly  reduced 
in  the  propane  flame  due  to  the  slower  flame-propagation  speed,  as  shown  in 
Fig.  8. 
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Table  1 


FUEL  AND  AIR  SETTINGS*  FOR  ACETYLENE  AND  PROPANE  FOR 
THE  PREMIXED  BURNER 


Inner  Air 
Nebulizer 
Outer  Air 
Inner  Fuel 
Outer  Fuel 
Argon  Sheath 


Acetylene 
9-14 
7  -  11 

28  -  40 

2.5  -  5 

4-9 
15 


Propane 
4-8 
3-7 
8-14 
0.2  -  0.5 

0.4  -  1 

15 


*A11  rates  are  in  liters/min. 


3.2  COMPUTER  CODES 


In  computer  codes  developed  to  date,  major  emphasis  has  been  placed  on  the 

determination  of  temperatures  as  represented  by  experimentally  obtained 

spectra.  The  calculation  of  the  theoretical  spectrum  is  an  important  part 

of  this  process  and  will  be  discussed  first.  To  complete  the  temperature- 
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determination  process,  a  nonlinear  Jacobian  least-squares  routine  has  been 
developed  for  fitting  the  observed  spectrum  as  a  function  of  temperature. 

The  third-order  susceptibility  which  governs  CARS  can  be  expressed  in 
abbreviated  form  by 


X  =  *n  +  *r  (ID 

where  y  is  the  nonresonant  contribution  and  y  Is  the  Raman  resonant 
n  r 

contribution  to  the  third-order  susceptibility. 
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(12) 


is  given  by 


X 


r 


-  2 

V,J 


6(v,j)  -  ir  (v,J) 


where  N  is  the  total  number  density  of  the  species  being  probed;  A(v,J)  is  the 
populational  difference  between  the  states  involved  in  the  Raman  transition; 
6(v,J)  is  the  detuning  factor  defined  by  S(v,J)  =  a>r  -  A  (w  being  the  fre¬ 
quency  of  the  Raman  transition  and  A  the  frequency  difference  between  the  pump, 
and  the  probe,  ii^) ;  T(v,J)  is  the  HWHM  of  the  Raman  transition;  and 
(do / dSJ) v  j  is  the  Raman  scattering  cross  section  defined  for  Q  branches  as 


\ 
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~ 2  +ihj  2" 

) 

'  V,J 
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x  7  o 

La  +  45  bJ  ^  _ 
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where  a  is  the  spherical  part  of  the  polarizability  and  bjy  is  the 
anisotropic  part  which  usually  has  a  weak  J  dependence. 
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The  CARS  power  is  given  by 

P3(u>)  =  J  P1(to1)dui1  |pi(w  -  A)P2(W;l  -  A)|X(A)|2dA  (13) 

In  the  normal  course  of  generating  a  CARS  spectrum,  the  third-order  suscepti- 

bility  is  generated  according  to  Eqs.  (11)  and  (12)  and  then  |xV  |  is 

convolved  over  the  pump  and  probe  laser  linewidths  as  dictated  by  Eq.  (13). 

Another  convolution  over  the  detector  slit,  however,  must  be  made  before 

the  calculated  spectrum  can  be  compared  to  the  experimentally  observed  one. 

Under  typical  experimental  conditions  when  probing  f^,  simplification  of 

Eq .  (11)  and  (12)  can  be  accomplished  without  the  introduction  of  large 

errors.  First,  the  nonresonant  third-order  susceptibility,  x  >  can  be 

n 

neglected,  except  for  the  cases  of  minor  species  concentrations.  Secondly, 
the  Raman  scattering  cross  section,  (dcr/dfi)  is  observed  to  be  a  slowly 

V  f  vl 

varying  function  of  J  level  and  can  be  treated  as  a  constant.  Thirdly,  the 
Raman  linewidths  vary  slowly  with  J  level  and  can  be  treated  as  constant 
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for  all  transitions.  Optimum  results  are  achieved  by  employing  a  linewidth 
of  0.1  cm  ^  for  the  spectral  calculations .  As  indicated  by  Eq.  (12),  the 
calculation  of  the  third-order  susceptibility  as  a  function  of  frequency 
involves  the  summation  over  all  Raman  resonances  in  the  spectrum.  For 
typical  Q-branch  spectra,  this  entails  the  summation  over  three  vibrational 
transitions,  each  containing  100  rotational  transitions.  Because  of  the 
narrow  Raman  linewidths,  the  frequency  increment  for  the  calculation  is 
kept  at  0.01  cm  \  covering  a  region  from  2265  to  2350  cm  ^ .  The  third- 
order  susceptibility  must  then  be  squared  and  convolved  over  the  appropriate 
experimentally  defined  slit  functions. 


(3)  2 

Figure  9  shows  the  |x  |  for  Nj  at  1700  K  before  convolution  over  the 
detector  halfwidth  and  laser  linewidths.  The  spectral  constants  employed 
for  this  calculation  are  given  in  Ref.  15.  The  present  program  has  the 


ability  to  convolve  over  any  of  four  different  functions  any  number  of 


times.  Figures  10  -  13  show  the  convolved  spectrum  utilizing  the  four 
available  functions  at  a  HWHM  of  1.67  cm  ^ .  The  formulae  for  the  individual 


functions  are: 


Gaussian 
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Figure  11.  A  Lorentzian  Slit  Convolution  of  ,  «  v  ;  at  1700 
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where  S  is  the  slope  of  the  rising  edge  of  the  trapezoid. 

The  experimental  setup  dictates  the  actual  form  which  Eq.  (13)  must  take  in 
order  to  simulate  the  obtained  spectrum.  In  the  laboratory  ruby  scheme,  the 
linewidth  of  the  pump  is  ~  0.03  cm  ^ ,  and  the  broadband  dye  has  a  linewidth 
of  '  200  cm  1.  The  CARS  signal  is  detected  by  an  JMA  which,  along  with  the 
monochromator,  contributes  to  a  1.2  c.m~^  halfwidth.  Under  these  experi¬ 
mental  conditions,  the  convolutions  over  the  pump  and  probe  lasers  can 
be  neglected  since  in  the  former  case  the  bandwidth  is  so  small  that  it 
has  little  effect  upon  the  spectrum;  and  in  the  latter  case  the  bandwidth 
is  so  large  that  it  is  in  constant  over  the  60  cm  ^  region  responsible  for 
CARS  production.  The  detector  slit  function,  however,  cannot  be  neglected 
and  must  be  convolved  over.  A  triangular  slit  function  is  primarily  employed 
for  this  convolution.  In  the  Nd:YAG  scheme  the  pump  has  a  bandwidth  of 
0.8  cm  ^  which  must  be  convolved  over;  the  dye  laser  has  a  bandwidth  of 
130  cm  ^  as  shown  in  Fig,  14  which,  as  in  the  case  of  the  ruby  scheme,  is 
not  convolved  over.  With  the  Nd’.YAG  system,  the  detector  halfwidth  increases 
to  2.0  cm  ^  due  to  the  decreased  dispersion  of  the  grating  in  this  frequency 
region.  According  to  Eq.  (13)  the  pump  bandwidth  must  be  convolved  over 
twice,  but  in  reality  a  single  Gaussian  convolution  was  found  to  be  adequate. 


A  nonlinear  least-squares  program  has  been  developed  to  fit  the  observed 
experimental  spectrum  as  a  function  of  temperature.  It  is  based  upon  the 
matrix  equation 

AP  =  (J+J)_1  J+  A<J> 
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Figure  14.  Broadband  Dye  Output  of  the  Nd : YAG-Pumped 
Dye  Laser. 


<1 


where  J  Is  the  Jacobian  matrix  whose  elements  are  defined  by  J.,  =  dl ./dP, 

th  th  +  ^  ^ 

(J^  being  the  i  intensity  and  the  k  parameter);  J  is  the  transpose 

of  the  J  matrix;  A<j>  is  the  matrix  with  elements  A«f> .  =  I,  ,  -■  I.  ,  (where 

t}i  i  iobs  .  ical 

^iobs  *s  t*ie  *  observed  experimental  intensity  and  l^ca^  the  ith  calculated 
intensity);  and  AP  is  the  correction  matrix  whose  elements  contain  the  cor¬ 
rections  to  be  made  upon  the  parameters  which  were  varied.  Three  parameters 
were  chosen  for  fitting  by  the  program — the  temperature,  the  relative  frequency, 
and  the  Raman  linewidths.  The  frequency  fit  was  incorporated  in  order  to 
negate  the  need  for  absolute  frequency  Calibration,  and  the  Raman  linewidth 
was  used  to  take  into  account  any  dependence  of  the  linewidth  upon  temperature. 


The  program  works  in  the  following  manner.  First,  the  initial  guesses  as  to 
the  parameters  to  be  varied  are  made  and  a  spectrum  calculated.  This  spectrum 
is  then  used  to  formulate  the  A$  matrix.  Secondly,  the  Jacobian  matrix  is 
created  by  a  numerical  evaluation  of  the  appropriate  derivatives.  The  trans¬ 
pose  of  the  Jacobian  and  subsequential  inverse  of  the  product  of  the  Jacobian 
with  its  transpose  is  formed.  The  correction  matrix  can  then  be  constructed 
according  to  Eq.  (18).  The  correction  to  each  parameter  can  then  be  made  and 
the  resulting  new  parameter  used  as  the  initial  guess  for  which  cycle  to  repeat. 
The  iterative  process  is  continued  until  the  variance  of  the  calculated  spectrum 
shows  no  further  improvement.  The  dispersion — an  indicator  for  the  fitting 
error  for  each  parameter — can  be  obtained  from 


wXJ  „2>\ 1/2 


(NOBS-NPAR) 


7 


(19) 


where  (J+J)^  is  the  it'1,  ith  diagonal  element  of  the  (J+J)  ^  matrix;  a ^  is 
the  variance  obtained  between  the  observed  and  calculated  spectra;  NOBS  is 
the  total  number  of  data  points  in  the  observed  spectrum;  and  NPAR  is  the 
total  number  of  parameters  varied.  The  iteration  usually  requires  three 
cycles  for  fitting  the  temperature.  (A  complete  listing  of  the  routine  which 
operates  on  a  CYBER  74  computer  can  be  found  in  the  Appendix.)  The  first  test 
of  t lit*  program  was  made  on  room-temperature  data  since,  due  to  the  overall 
lack  of  structure  in  the  spectra,  it  is  probably  the  worst  possible  case  to 
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fit.  These  fits  were  routinely  within  15  K  of  room  temperature.  One  such 
fit  is  shown  in  Fig.  15.  Figure  16  depicts  a  fit  that  was  made  to  single¬ 
shot  data  obtained  from  the  premixed  propane  burner  by  the  ruby  system.  The 
experimental  data  obtained  from  the  same  burner  by  the  Nd:YAG  CARS  system  are 
shown  in  Fig.  17. 

3.3  EFFECT  OF  LINEWIDTHS  UPON  TEMPERATURE  MEASUREMENTS 

As  outlined  above,  the  temperature  of  a  combustion  medium  can  be  estimated 
by  obtaining  the  best  fit  of  a  calculated  spectrum  to  an  experimentally 
determined  one.  To  achieve  this  fit,  temperature  and  linewidth  can  be  used 
as  adjustable  parameters  if  the  J  dependence  upon  the  linewidth  in  the 
Q-branch  is  neglected.  This  approach  has  yielded  excellent  results  for 
laminar  laboratory  flames  where  the  temperature  has  been  determined  by  other 
methods  such  as  the  line-reversal  technique. However,  it  would  be  advan¬ 
tageous  to  gain  a  better  understanding  of  the  linewidth  problem  to  improve 
the  reduction  of  CARS  data  to  temperatures  and  species  concentrations.  The 
latter  case  presents  a  much  more  difficult  problem  and  will  be  discussed 
later. 

The  primary  mechanisms  for  spectral  line  broadening  in  a  combustion  medium 

are  molecular  collisions  and  the  Doppler  effect.  The  linewidttys  depend  upon 

temperature,  pressure,  and  composition  of  the  medium  as  well  as  the  J  value 

of  the  line.  Collisional  effects,  as  demonstrated  for  the  case  of  hydrogen, 

may  cause  narrowing  or  broadening  of  the  spectral  line.^  Fortunately, 

collisional  narrowing  can  be  neglected  for  most  species.  However,  it  does 

present  a  special  problem  for  the  analysis  of  hydrogen  data,  and  additional 

work  is  required  for  proper  analysis  of  data  in  combustion  environments. 

As  an  approximation,  the  collisional  linewidth  can  be  expressed  as  a  function 

18 

of  the  optical  collisional  frequency,  which  yields,  for  the  half-intensity 
breadth, 

T  =  ~  (20) 
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Figure  15.  Room-Temperature  Fit  of  OMA  Data. 
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Figure  16.  Temperature  Fit  of  a  Single-Shot  Spectrum  Obtained  with  the  Ruby-Based  CARS 
System  on  the  Propane-Fueled  Burner.  The  solid  line  is  the  calculated  fit, 
and  the  dots  are  the  experimentally  observed  points. 


where  is  the  time  between  collisions.  From  kinetic  theory  one  obtains, 
for  a  gas  composed  of  a  single  species. 


r 


2 


P 


(21) 


3 

where  p  is  the  optical  collisional  diameter  which  may  be  10  to  10  times 

greater  than  the  value  obtained  from  kinetic-theory  experiments.  For  N 

-3  ^ 

this  yields  an  estimated  half-intensity  breadth  between  2  *  10  and 

2  x  10  1  cm  1  at  atmospheric  pressure  and  room  temperature.  If  one  assumes 

that  the  collisional  diameter  is  independent  of  the  velocity  of  the  molecules, 

-1/2 

then  T  should  depend  upon  T  .  This  assumption  appears  to  be  in  reasonable 

19 

agreement  for  linewidth  measurements  using  CARS.  A  review  of  recent 

experiments  and  theoretical  developments  in  the  area  of  Raman  line  shapes  is 

20 

given  by  Srivastava  and  Zoidi.  The  theory  was  compared  with  experimental 

21 

CARS  data  for  the  N2  Q-branch  by  Rahn,  Owyoung,  Coltrin,  and  Koszykowski. 
Although  these  data  indicate  that  the  J  dependence  is  weak,  a  sufficient 
change  in  V  occurs  to  warrant  the  inclusion  of  this  effect  in  calculated 
spectra.  Additional  work  is  required  to  investigate  the  effect  of  neglecting 
the  J  dependence  when  temperatures  are  determined  by  fitting  measured  Q-branch 
spectra  to  calculated  spectra. 


As  the  temperature  of  the  combustion  gas  increases  for  a  given  pressure,  the 
effect  of  Doppler  broadening  becomes  more  important.  Thus,  the  linewidth 
becomes  a  combination  of  Doppler  and  collisional  broadening.  If  collisional 
narrowing  is  neglected,  it  can  be  shown  that  the  anti-Stokes  susceptibility 
becomes ^ 
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where 


NA  c 

K.  T  =  — V 

V*J  3hW4 
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and  AWp  is  the  Doppler  width  for  forward  Raman  scattering  given  by 


Au)^  =  (Dj[8(ln2)kT/mc2]^2 


(23) 


If  one  lets 
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2(ln2) 
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the  equation  for  x 
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where  z  =  [  (cd  - 
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If  T  ->  0,  it  can  be  shown  that  in  this  Doppler  limit,  the  expression  reduces 
to 
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where 


dt 


1/2 

(Dawson's  integral)  is  tabulated  in  Ref.  22  and  z  =  2(ln2)  (w  -  w  )/Auip. 
The  effect  of  the  real  part  is  to  increase  the  halfwidth  of  the  Doppler  line 
and  to  increase  the  wings.  At  first,  the  real  part  was  estimated  to  be 
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small.  With  this  assumption  the  width  of  |x'  ; (w) |  in  the  Doppler  limit 

becomes  Later  it  was  shown  that  the  real  part  does  contribute 

significantly  and  that  the  halfwidth  is  -  1.22  Aw  for  a  Doppler-broadened 
23  D 

CARS  line. 

If  Au>d  -*■  0,  the  Gaussian  function  in  Eq.  (22)  becomes  a  6-function,  and  one 

obtains  the  ordinary  collisional  broadened  line.  However,  in  high-temperature 

regions  of  a  combustion  medium,  Doppler  broadening  can  become  significant. 

_  o 

For  example,  at  2000  K  the  Doppler  width  for  becomes  2.4  x  10  which  is 
comparable  to  the  collisional  width  at  atmospheric  pressure.  The  J  dependence 
for  the  total  linewidth  would  be  reduced  under  these  conditions.  When  it  is 
necessary  to  include  Doppler  broadening,  |x'  ' (w) |  can  be  estimated  by  using 
Eq.  (22),  if  the  nonresonant  background  is  not  significant.  When  the  J 
dependence  is  known,  then  the  temperature  can  be  determined  by  fitting  cal¬ 
culated  data  to  experimental  data  using  temperature  and  T  as  adjustable 
parameters. 

3.4  LASER-INDUCED  GAS-BREAKDOWN-THRESHOLD  MEASUREMENT 

As  with  many  physical  processes  initiated  by  laser  radiation,  the  intensity 
of  the  radiation  generated  in  the  CARS  process  by  mixing  two  laser  beams 
depends  critically  upon  the  intensity  of  the  incident  laser  beams  in  an 
interaction  medium.  In  general,  this  means  that  the  more  intense  the  incident 
beams,  the  more  intense  the  generated  CARS  radiation.  There  is,  however,  a 
practical  .limit  beyond  which  an  increase  in  incident  intensity  will  not 
result  in  increased  generation  of  CARS  radiation.  One  of  the  processes  which 
impose  this  limit  upon  the  incident  laser  intensity  is  the  laser-induced  gas- 
breakdown  phenomenon.  Since  the  gas  breakdown  severely  disturbs  the  optical 
properties  of  the  medium  through  which  the  beams  propagate,  its  occurrence 
cannot  be  tolerated  in  optical-measurement  processes  such  as  CARS. 

In  order  to  determine  the  upper  limit  in  the  laser  intensity  (or  power)  which 
can  be  safely  employed  in  flame  diagnostics  without  causing  gas  breakdown, 
an  experiment  was  carried  out  to  measure  the  breakdown  threshold  intensity 
for  the  flat  flame  of  the  Battelle  burner. 
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The  experimental  set-up  is  shown  in  Fig.  18.  A  single-mode  ruby  laser  beam 
is  generated  in  an  oscillator  cavity  and  then  passed  through  two  amplifiers. 

A  portion  (~  0.1%)  of  the  beam  is  detected  by  a  PIN  photodiode  and  displayed 
on  a  Tektronix  R7912  transient  digitizer.  The  majority  of  the  ruby  beam 
continues  on  through  a  collection  of  colored-glass  and  neutral-density  filters 
and  enters  a  glass  cross  through  a  3-in.-diam.  quartz  window.  Inside  the 
cross  the  ruby  energy  is  focused  by  a  55-mm-focal-length  quartz  lens.  The 
focal  point  of  this  lens  is  near  the  vertical  centerline  of  a  2. 5-in.-diam. 
low-pressure  flat-flame  burner  developed  by  Battelle  Laboratories.  The  laser 
beam  exits  the  cross  through  another  3-in.  quartz  window  and  is  detected 
through  a  Schott  KG1  filter  by  means  of  a  thermopile  detector  head  which  is 
connected  to  a  Scientech  362  power /energy  meter. 

Two  steps  were  taken  to  minimize  fluctuations  in  the  energy  meter  readings 
caused  by  environmental  conditions  and  background  radiation.  The  effects  of 
temperature  gradients  in  the  surrounding  air  upon  the  thermopile  were  reduced 
by  use  of  an  insulated  enclosure  which  allowed  a  l-in.-diam.  access  to  the 
detector  head.  This  access  port  was  covered  with  a  Schott  KG1  filter  which 
was  used  to  absorb  the  IR  radiation  produced  in  the  flame.  With  this  filter¬ 
ing  background  energy  affecting  the  detector  was  independent  of  the  observed 
position  within  the  flame. 

The  experimental  procedure  used  to  arrive  at  the  laser  energy  required  for 
media  breakdown  consisted  of  monitoring  the  laser  pulse  shape  on  the  transient 
digitizer,  recording  the  energy-meter  reading,  and  adjusting  the  laser  power 
to  the  required  levels.  The  PIN  diode/digitizer  combination  supplied  infor¬ 
mation  on  the  pulsewidth  and  the  relative  amplitude  of  the  laser  pulse. 

The  optical  breakdown  of  the  gaseous  media  induced  by  the  laser  radiation 
manifests  itself  in  three  ways,  the  first  two  being  the  generation  of  a 
visible  discharge  and  an  audible  "snap"  in  the  focal  volume.  The  third 
phenomenon,  which  was  most  heavily  relied  upon,  was  the  pronounced  reduction 
in  the  laser  energy  transmitted  through  the  affected  media,  the  media  break¬ 
down  reducing  by  a  factor  of  two  or  more  the  energy  throughput. 
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Figure  18.  Schematic  Diagram  of  the  Experimental  System 
for  Breakdown-Threshold  Measurement. 
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Adjustment  of  the  amount  of  laser  energy  being  focused  in  the  flame  cell  was 
accomplished  in  two  ways.  Large-scale  power  changes  were  accomplished  by 
insertion  of  various  combinations  of  neutral-density  filters  into  the  beam. 
Fine  adjustments  in  intensity  were  made  by  changing  the  ruby-amplifier 
flash-lamp  voltage.  Thus,  by  increasing  the  laser- pulse  energy  to  the 
breakdown  limit  and  measuring  the  energy  transmitted  in  a  "clean"  Gaussian 
beam,  the  point  of  optical  breakdown  was  determined. 

To  arrive  at  the  power  density  required  to  cause  breakdown  in  the  media,  the 
following  expression  was  employed: 


d  T  T  y  (1.27  A  f/D)2  (26^ 

pw  c  4 

where  P  =  power  density;  E,  _  =  detected  laser  energy;  r  =  laser  pulse- 
d  det  pw 

width  (FWHM);  Tc  =  combined  transmission  of  exit  window  and  KG1  filter;  A  = 

ruby-laser  wavelength;  f  =  focal  length  of  focusing  lens;  and  D  =  diameter  of 

laser  beam  at  focusing  lens.  The  values  of  the  parameters  used  in  these 

measurements  were:  t  =13  ns,  T  =  62.6%,  A  =  6943  A,  f  =  53  mm,  and  D  = 

pw  c 

3.5  mm. 


Figure  19  shows  the  variation  of  breakdown  power  densities,  as  determined 
through  the  above  procedure,  exhibited  for  five  sets  of  conditions  existing 
at  the  focal  volume.  Curves  1  and  2  show  the  breakdown  characteristics  at  1 
atm  of  pressure,  while  3  and  4  show  the  effects  observed  at  130  Torr.  Data 
Point  5  indicates  the  trend  of  the  breakdown  threshold  in  a  fuel-rich  (dirty) 
flame. 

One  general  conclusion  that  is  supported  by  the  data  in  Fig.  19  is  that  the 
power  density  required  for  media  breakdown  is  inversely  proportional  to  the 
number  density  of  the  molecules  available  in  the  focal  volume  and  that  this 
number  density  can  be  varied  by  changing  either  the  pressure  or  temperature. 

A  comparison  of  Curves  2  and  4  lends  credence  to  this  conclusion  and  reflects 
known  characteristics  of  flame  composition  and  geometry.  Specifically,  the 
leading  edges  of  these  two  curves  tend  to  describe  a  mixing  zone  for  the 
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Figure  19.  Laser-Induced  Breakdown-Threshold  Profile  of  the  Battelle-Burner  Flame. 


reiativcj.^  ~jo1  combustion  gases.  Also,  the  trailing  edge  of  these  curves 
describes  the  effect  of  pressure  upon  a  flame  operating  in  a  lower-pressure 
environment  causing  a  lengthening  of  the  flame  reaction  zone.  To  our  relief, 
the  reduction  in  breakdown  threshold  due  to  high  temperature  in  flames  which 
had  been  anticipated  did  not  materialize  in  our  experiments. 

During  the  course  of  this  experiment,  one  heretofore  unnoticed  physical 
characteristic  of  the  Battelle  burner  was  observed.  The  surface  of  the 
burner  consists  of  a  2. 5-in.-diam.  wire  mesh  which,  when  heated,  buckles 
slightly.  When  breakdown  measurements  were  attempted  at  altitudes  less  than 
1  nun  above  the  burner,  the  wire  screen  intercepted  the  laser  radiation, 
resulting  in  undersirable  scattering  and  absorption.  This  characteristic,  in 
effect,  limits  the  usable  measurement  zone  of  this  burner  to  a  minimum  of 
1  mm. 


3.5  POSSIBILITY  OF  USING  THE  CO  RAMAN  BANDS  FOR  REAL-TIME  TEMPERATURE 
MEASUREMENT  IN  FLAMES 

Efforts  on  developing  a  temperature-measurement  technique  for  application  to 

flame  diagnostics  using  the  CARS  process  have  thus  far  been  limited  to  the 

analysis  of  the  band  spectrum  in  one  manner  or  another.  This  analysis  may 

involve  either  evaluation  of  the  slope  of  the  linear  portion  of  the  Q-branch 
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spectrum  as  outlined  previously,  or  contour  matching  of  an  experimental 
spectrum  with  those  generated  by  a  computer.  An  alternate  approach  to  obtain¬ 
ing  temperatures  from  combustive  environments  is  to  use  the  vibrational 
structure  of  the  CARS  spectrum  of  the  CO^  molecule  present  in  flames.  This 
method  is  based  upon  the  comparison  of  integrated  intensities  for  one  of  the 
Raman-active  fundamentals  and  its  associated  hot  bands.  This  method,  if  it 
proves  to  be  practical,  will  have  a  number  of  advantages  over  use  of  the  N^ 
spectrum.  First,  because  this  method  is  insensitive  to  fluctuations  in  the 
CO^  concentration  in  flames,  it  would  permit  real-time  as  well  as  time- 
averaged  temperature  measurement.  Secondly,  the  method  would  require  a 
relatively  simple  data-processing  system  and  instrumentation  for  operation. 

The  disadvantages  are  that  the  temperature  measured  is  vibrational  temperature 
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of  a  combustion  product  and  thus  may  be  reaction-dependent  and  that  the 
relative  concentration  of  CO 2  available  in  flames  typically  at  -  10%  is  lower 
than  that  of  N2,  and  thus  may  be  difficult  to  detect. 


Although  the  CO^  molecule  has  only  one  Raman-active  fundamental  band,  due  to 
Fermi  resonance*  ,  there  are  in  reality  two  strong,  fundamental-like  bands 
at  1388.2  cm  ^  and  1285.4  cm  ^  which  correspond,  respectively,  to  (10%  -  00%) 
and  (02%  -  00%)  transitions.  Although  the  assignment  used  here  is  the 
traditional  assignment  of  these  bands,  it  has  been  shown  in  recent  years 
that  this  traditional  assignment  should,  in  fact,  be  reversed  (see  Refs.  26 
and  27).  Each  of  these  bands  has  its  own  hot-band  structure  which  becomes 
more  prominent  as  the  gas  temperature  is  increased.  The  pertinent 
energy  levels  and  the  Raman  bands  of  interest  are  diagramatically  shown  in 
Fig.  20.  For  diagnostic  purposes,  either  of  the  two  bands  can  be  used. 
However,  the  one  labeled  1388  cm  ^  group  is  about  twice  as  strong  as  the 
other,  and  naturally  it  is  a  better  choice.  Also,  the  1388  cm  ^  group  does 
not  suffer  from  interference  of  the  other  blending  bands  as  does  the  1285  cm 
group.  For  example,  the  Raman  transition  (04%  -  02%)  at  1263.0  cm  * 
(although  not  shown  in  the  diagram)  is  very  near  the  (03%  -  01%)  band  at 
1265.1  cm  %  Since  the  relative  intensity  of  these  lines  is  a  function  of 
the  population-density  difference  between  levels  involved  in  the  CARS  process 
and  since  the  population  density  is,  in  turn,  a  function  of  temperature, 
measurement  of  the  relative  intensity  of  two  or  more  lines  can  be  used  for 
determining  the  gas  temperature. 


The  integrated  power  of  the  first  hot  band  at  1409.5  cm  *  relative  to  that 
of  the  fundamental  band  at  1388.2  cm  can  be  expressed  by: 


PINT,1 

PINT,0 


FX(T) 

Fq(T) 


reff,0 

Feff,l 


(27) 
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Figure  20.  Pertinent  Energy  Levels  and  Raman  Bands  of  Interest  to 
Combustion  Diagnostics.  Energy  levels  (in  cm-^  units) 
were  obtained  from  Ref.  28. 


where  is  the  integrated  CARS  power  of  a  band,  is  the  number-density 

difference  between  the  lower  and  upper  levels  involved  in  the  CARS  process, 
d©/df2  is  the  differential  spontaneous  Raman  cross  section,  and  F(T)  and  F  ^ 
are  the  previously  defined  parameters^  representing  the  effect  of  rota¬ 
tional  distribution  and  the  interference  among  various  rotational  lines, 
respectively.  In  the  following  paragraphs,  the  factors  which  appear  in 
Eq.(27)  will  be  examined  individually. 


In  general,  the  population  density  of  a  vibrational  level  of  a  linear  triatomic 
molecule  (ignoring  the  vibrational  angular  momentum  and  Fermi  resonance  for  the 
moment)  is  given  by: 


»‘V  W  ■ 


exp 


-E(v1*v2,V3) 

kT 


(28) 


where  v^.v^  an<*  V3  are  t*ie  three  vibrational  quantum  numbers  of  the  C02 
molecule;  ECv^.Vj.v^)  is  the  energy  of  the  level  (v^.v^v^);  and  Qv  is 
the  vibrational  partition  function  given  by: 


VT>-  2 1  2(v2  +  x)  exp 

V1,V2’V3 

where  the  temperature  dependence  of  the  partition  function  is  explicitly 
indicated. 


-E(v1,v2,v3) 


(29) 


Vibrational  levels  with  v2  >  2  are  split  according  to  their  vibrational 
angular  momentum  which  takes  on  a  value  ~  v2>  v2~2,...l  or  0.  The 
statistical  weight  of  these  sublevels  is  a  function  of  and  is  tabulated 
below  for  the  first  three  excited  levels  of  v2. 
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A  real  CARS  spectrum  ot  the  1366  cm  A  group  is  shown  in  Figv..  ?t.  This  s, 
trum  having  a  spectral  resolution  of  2  cm""1  was  recorded  with  an  OMA  from  the 
flat  flame  of  the  Battelle  burner  operated  with  the  N2  buffer  gas  replaced 


by  CO 2  gas.  Note  that  the  lines  at  1388.2  and  1409.5  cm"1  were  attenuated 
by  a  factor  of  10  (ND  1)  and  4  (ND  0.6),  respectively.  Even  at  room  tempera¬ 
ture,  the  first  hot  band  (11*0  -  01*0)  can  be  easily  observed  since  about 
7%  of  the  CO 2  population  occupies  the  (01*0)  state. 


The  ratio  of  the  number-density  difference  of  the  first  excited  state  (01*0) 
to  that  of  the  ground  state  (00°0)  is  given  by: 


Ndl  _  N(011Q)  -  N(11*0) 

NJn  N(00°0)  -  N(10°0) 
au 


l 

N(01  0) 

_  1  -  exp 

E(11*0)  ~  E(01*0)\ 

V'  .  W . j, 

— 

N(00°0) 

1  -  exp 

(  E(10°0)  -  E(00°0)\ 

V  kT  / 

»  ,  2  exp  F-  960. 0/t]  (3f)> 

N(00°0)  L  J 


where  N(00°0)  is  the  ground-state  number  density  and  the  others  are  similarly 
defined.  The  second  equality  is  justified  in  view  of  the  nearly  equal  energy 
differences  in  the  exponents.  This  approximation  entails  less  than  0.3%  error 
at  2000 °K ,  although  at  that  temperature  the  upper- level  population  is  more  than 
30%  that  of  the  lower  level.  Similarly  for  the  second  excited  state,  the  ratio 
can  be  written  aa: 


Nd2  N(Q22Q)  -  N(1220)  N(0220) 

N~  s  N(00°0)  -  N(10°0)  “  N(00°0) 


2  exp 


-1920. 5/T 


(33) 


These  ratios  are  plotted  in  Figure  22  as  a  function  of  temperature.  One 
measure  of  the  sensitivity  of  the  temperature  measurement  is  given  by  the 
derivative  of  these  ratios  with  respect  to  temperature.  However,  the  accuracy 
with  which  the  population  ratio  can  be  determined  in  a  CARS  experiment  is 
usually  given  in  terms  of  a  percentage  error.  Therefore,  it  is  more  meaning¬ 
ful  in  practice  to  consider  the  percentage  change  in  the  population  density 
ratio  with  temperature  than  the  simple  derivative  in  this  regard.  This 
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Figure  21.  CARS  Spectrum  of  C0„  in  a  Flame  in  the  Virinit 
of  1400  cm-1.  1 
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Ratio  of  the  Population  Densities  of  (01-*-0)  and  (02^0)  States  of 
the  CO2  Molecule  to  that  of  the  Ground  State  (OOOO).  Preliminary 
experimental  results  are  also  shown. 


percentage  change  is  given  by 


d_ 

dT 


960. 0/T2 


for  the  first  hot  band  and 


d_ 

dT 


1920. 5/T2 


C2) 


(3d) 


for  the  second  hot  band.  This  indicates  that  the  sensitivity  of  the  second 
hot  band  is  greater  than  that  of  the  first.  In  fact,  the  sensitivity  is 
directly  proportional  to  the  energy  spacing  of  the  levels  being  compared. 

Of  course,  this  criterion  of  the  sensitivity  is  valid  only  as  long  as  the 
population  densities  can  be  determined  with  the  same  accuracy  for  both 
levels . 


Equation  (33),  for  instance,  can  be  used  to  estimate  the  accuracy  of  the 

temperature  measurement  made  using  the  second  hot  band  and  the  fundamental 
band.  For  this  estimation  it  is  more  convenient  to  rewrite  the  equation 
as  follows: 


dT 

T 


(34) 


The  above  equation  states  that  the  estimate  of  the  error  in  the  temperature 
measurement  is  T/1920.5  times  the  relative  error  in  the  ratio  of  the  number 
densities.  Thus,  a  5%  error  in  the  determination  of  the  ratio  of  the  number 
densities  would  result  in  a  4.7%  and  2.6%  error  at  1800  K  and  1000  K,  respec¬ 
tively.  In  comparison,  the  use  of  the  first  hot  band  in  conjunction  with 
the  fundamental  would  yield  9.4%  and  5.2%,  respectively, for  the  same  tempera¬ 
tures  . 


Function  F(T) ,  which  accounts  for  the  effect  of  changing  rotational  distri¬ 
bution  with  temperature,  has  been  calculated  for  the  first  two  excited 
1  2 

states  (01  0  and  02  0)  as  well  as  the  ground  state,  and  the  results  are 
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tabulated  in  Table  2.  Neglecting  to  include  this  factor  in  the  evaluation  of 

the  number  density  from  CARS  intensity  data  would  introduce  an  error  of  up  to 

1  ? 

3%  and  6%  for  the  (01  0)  and  (02  0)  levels,  respectively. 

Information  on  the  effective  linewidrh  (T  , , )  and  the  differential  Raman  cross 

eff 

section  for  the  Raman  bands  under  consideration  is  not  available  at  this  time. 

Pressure-broadened  Raman  linewidth  data  are  needed  to  calculate  I'  .  However, 

ef  f 

because  of  the  small  difference  in  the  rotational  constants  of  the  states 
involved  in  the  CARS  process,  the  Q-branch  lines  are  not  resolved  at  all,  and 
thus  it  would  be  impossible  to  measure  linewidth.  It  would  be  highly  desirable 
to  measure  the  Raman  cross  section  of  the  hot  bands  relative  to  that  of  the 
fundamental  band. 

Preliminary  results  of  experiments  designed  to  test  the  feasibility  of  this 

idea  were  also  presented  in  Fig.  22.  The  CO^  gas  was  heated  in  an  oven  and 

the  temperature  monitored  with  a  thermocouple.  The  integrated  CARS  power  was 

measured  using  a  Stokes  beam  having  a  linewidth  of  ~  2  cm  1 .  This  linewidth 

is  sufficiently  wide  to  excite  all  the  rotational  lines  that  are  thermally 

populated  since  calculation  shows  that  the  widths  of  the  Q-branches  are  in 

fact  0.10,  0.75,  and  0.88  cm"1  for  the  (10°0  -  00°0),  (1110  -  0110)  and 
2  2 

(12  0  -  02  0)  bands,  respectively.  In  making  these  measurements,  the  Stokes 

beam  was  tuned  from  band  to  band  at  each  temperature  setting.  A  least-squares 

fit  of  the  experimental  data  to  the  theoretical  curve  yielded  a  value  of  0.31 

for  the  ratio  of  the  factor  (do/dfl)//rT7  for  the  first  hot  band  to  that  of 

eff 

the  fundamental  and  similarly  0.49  for  the  second  hot  band.  Although  the 
present  experi  nt  is  not  being  carried  out  in  real  time,  the  extension  of 
this  scheme  to  real-time  measurement  should  not  be  very  difficult. 
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Table  2 


CALCULATED  VALUES  OF  F(T)  FOR  C02  NORMALIZED  TO  THE  VALUE  FOR  T 


T 

Fo(T) 

F1(T) 

F2(T) 

(K) 

Fo(300) 

FiOOO) 

F2(300) 

300 

1.000 

1.000 

1.000 

450 

0.819 

0.815 

0.815 

600 

0.711 

0.706 

0.707 

750 

0.635 

0.630 

0.630 

900 

0.579 

0.575 

0.575 

1050 

0.536 

0.529 

0.529 

1200 

0.503 

0.493 

0.494 

1350 

0.479 

0.464 

0.464 

1500 

0.462 

0.441 

0.441 

1650 

0.449 

0.424 

0.419 

1800 

0.440 

0.411 

0.398 

1950 

0.432 

0.403 

0.381 

2100 

0.428 

0.397 

0.366 

300  K 


3. b  COMPARISONS  OF  SINGLE-SHOT  THERMOMETRY  OF  CARS  WITH  OTHER  OPTICAL 
THERMOMETRIC  TECHNIQUES 


Two-Line  Fluorescence 


A  fluorescence  system  was  developed  in  order  to  implement  the  two-line  fluo- 
29 

rescence  techniques  for  measuring  flame  temperature.  The  development  of 
this  system  was  motivated  by  the  desire  for  a  means  of  determining  temperature 
in  flames,  independent  of  the  CARS  technique.  An  independent  measurement  of 
flame  temperature  would  constitute  a  firm  basis  for  evaluating  the  validity 
and  accuracy  of  the  CARS  process  as  a  thermometric  tool  for  flames. 

30 

Winefordner  has  shown  that  the  temperature  in  a  flame  can  be  measured  by 
means  of  a  two-line  fluorescence  method.  The  advantage  which  this  method  has 
over  other  spectroscopic  methods  that  have  been  used  is  spatial  resolution, 
i.e.,  "local"  temperatures  in  various  parts  of  the  flame  can  be  measured 
directly  rather  than  inferred  from  average  temperature  measurements  across 
the  flame.  The  method  involves  seeding  the  flame  with  atoms  whose  energy- 
level  structure  is  similar  to  that  shown  in  Fig.  23  (for  example.  In,  Ga,  or 
Tl).  Level  2  is  connected  via  electric  dipole  transition  to  both  Levels  0 
and  1.  It  is  assumed  that  thermal  equilibrium  exists  and  the  ratio  of  the 
populations  of  Levels  1  and  0  is  given  by 


N 

N 


1 

0 


=  —  exp(-E  /kT) 
g0  1 


(35) 


where  g^  and  are  the  statistical  weights  of  Levels  1  and  0,  respectively; 

E^  is  the  energy  difference  between  Levels  1  and  0;  k  is  the  Boltzmann  con¬ 
stant;  and  T  is  the  local  temperature  of  the  flame.  E^  is  generally  known. 
Therefore,  determination  of  the  ratio  N^/N^  would  yield  a  value  for  T.  A 
simple  analysis  shows  that  N^/N^  may  be  measured  in  the  following  way:  The 
transition  0  +  2  at  wavelength  is  excited  by  a  broadband  light  source  of 

n  —1. 

known  energy  density  (erg  cm--'  Hz  )  centered  at  an<*  F^  fluorescent 

photons/sec  are  observed  at  wavelength  Then  under  steady-state  conditions 
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Figure  23.  Level  Diagram  for  Two-Line  Fluorescence. 
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CN2  A21 


(37) 


where  C  is  a  constant  which  depends  upon  a  geometrical  factor,  and  Z^  is  the 
total  non-radiative  quenching  rate  of  Level  2.  The  stimulated -emission  term 
(Pq2®20N2^  w*1*c*1  is  very  small  has  been  neglected.  Also  the  transition  1  ■+  2 
is  excited  by  a  broadband  light  source  of  energy  density  centered  at 
wavelength  and  F^q  fluorescent  photons  are  observed  at  wavelength 

Under  steady-state  conditions  one  may  write  two  equations  which  are  similar 
to  Eqs.  (36)  and  (37).  A  simple  manipulation  of  these  four  equations,  along 
with  the  use  of  the  well-known  relationship  between  Einstein  A  and  B  coef¬ 
ficients,  yields 


F20  h 
F21  g0 


P02 

P12 


(38) 


Note  that  the  ratio  N^/N^  is  independent  of  the  quenching  rate  Z^,  even  under 
nonsaturating  conditions.  Using  the  measured  values  of  Pqj*  pj2*  F20’  ant* 

F21’  t*le  temPerature  T  may  be  determined  using  Eqs.  (35)  and  (38).  Winefordner 
has  measured  the  temperature  profiles  of  several  kinds  of  flames  using  this 
method. 


The  experimental  system  is  shown  schematically  in  Fig.  24.  The  broadband 
light  output  from  an  Eimac  lamp  is  chopped  (~  145  Hz),  line-selected  for 
either  410.2  or  451.1  nm  by  means  of  a  pair  of  interference  filters,  and 
focused  (LI)  into  an  argon-shielded,  premixed  acetylene  flame  which  is  seeded 
with  indium  atoms  (in  reality,  InCl^  dissolved  in  distilled  water  is  nebulized 
into  the  flame).  In  the  two- line  fluorescence  method,  the  fluorescent  signal 
is  measured  at  one  wavelength — say  410.2  nm — while  the  fluorescing  medium  is 
excited  at  the  other  wavelenght — 451.1  nm — and  vice  versa.  The  nonresonant 
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Figure  24  .  Schematic  Diagram  of  a  Two-Line-Fluorescence  System. 


fluorescence  signal  from  the  excited  indium  atoms  is  collected  (L2),  analyzed 
(monochromator),  and  synchronously  detected  (photomultiplier  plus  lock-in 
amplifier),  and  recorded.  At  the  same  time,  the  dc  component  of  the  photo¬ 
multiplier  current  is  monitored  and  recorded  in  order  to  ascertain  the  con¬ 
stancy  of  the  level  of  seeding  throughout  a  measurement . 

In  making  the  two-line  fluorescence  measurement,  it  is  necessary  that  the 
relative  transfer  function  (defined  here  as  the  ratio  of  the  photomultiplier 
output  to  the  incident-light  intensity)  be  known.  This  was  accomplished 
using  an  NBS-calibrated  standard  tungsten  filament  lamp.  Figure  25  gives  an 
example  of  such  a  transfer  function  for  the  detection  system  used  in  the 
experiment  which  consists  of  a  Spex  1702  monochromator  having  a  1200-gr/mm 
grating  and  an  EMI  9816  photomultiplier.  The  values  of  the  transfer  function 
for  the  two  fluorescence  lines  i.e.,  410.2  and  451.1  nm,  were  calculated  to 
be  2.4169  and  2.9134,  respectively,  and  their  ratio  was  0.8296.  A  spectral 
scan  of  the  transfer  function  with  a  different  standard  lamp  yielded  a  similar 
result  for  the  relative  transfer  function  at  the  two  wavelengths. 

In  order  to  determine  the  optimum  seeding  level  for  the  measurement,  a  set  of 
stock  solutions  having  various  concentrations  of  InCl^  was  prepared;  the 
fluorescent  intensity  was  measured  for  the  two  wavelengths  as  a  function  of 
InCl^  concentration.  The  results  which  are  shown  in  Fig.  26  indicate  that 
the  fluorescence  is  a  linear  function  of  the  InCl^  concentration  up  to 
200  mg/1.  At  higher  concentrations,  the  self-absorption  of  the  medium  can  no 
longer  be  neglected,  i.e.,  the  medium  cannot  be  treated  as  though  it  were 
optically  "thin."  Based  upon  these  observations,  an  optimum  concentration  of 
150  mg/1  was  used  for  all  experiments. 

Data  obtained  from  the  two-line  atomic  fluorescence  system  are  presented  in 
Fig.  27.  This  profile  was  obtained  by  stepping  the  approximate  6-mm-diam. 
focused  pump  beam  through  the  10-mm-diam.  indium-seeded  flame  at  a  height  of 
4  cm  above  the  burner.  The  profile  indicates  that  the  temperature  is  con¬ 
stant  over  a  large  section  of  the  burner.  The  measurement  resulted  in  a 
temperature  of  2650  +  50  K.  Several  difficulties  were  encountered  in  per¬ 
forming  the  two-line  fluorescence  experiments.  Of  most  concern  was  the 
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Relative  Transfer  Function  of  the  Detection  System 
(Monochromator  and  Photomultiplier). 


relatively  low  signal-to-noise  ratio  (3:1)  obtainable.  This  required  the  use 
of  long  lock- in-amplifier  filter  time  constants  (30  sec),  which  resulted  in 
slow  system  response. 

In  order  to  compare  the  results  to  CARS  measurements,  a  similar  profile  of 
the  acetylene-fueled  flame  was  obtained  with  the  Nd:YAG-based  CARS  system. 

This  profile  (shown  in  Fig.  28),  like  the  two-line  fluorescence  results, 
indicated  that  the  temperature  was  essentially  constant  over  a  large  region 
of  the  flame  and  dropped  sharply  at  the  edges.  The  average  temperature  from 
these  sets  of  measurements  was  2600  K  +  50  K.  This  was  generally  in  good 
agreement  with  the  fluorescence  measurement,  considering  the  larger  uncer¬ 
tainty  of  those  data.  Figure  29  shows  a  computer  fit  which  was  made  to  the 
data,  indicating  that  the  program  and  input  parameters  were  indeed  replicating 
the  true  CARS  spectra. 

Sodium-Line  Reversal 


A  system  for  performing  temperature  determinations  by  the  sodium  line-reversal 
31 

technique  was  assembled,  as  shown  in  Fig.  30.  This  technique  utilizes  the 
emission  characteristics  of  (in  this  case)  the  sodium  molecules  in  ordinary 
tap  water  nebulized  into  the  flame  to  determine  a  blackbody  radiance  function 
from  which  the  flame  temperature  is  extracted.  An  average  temperature  of 
2560  K  was  measured  for  the  6-mm-diam.  sample  volume  through  the  acetylene- 
fueled  flame  at  the  test  height  of  4  cm.  This  is  in  excellent  agreement  with 
the  CARS  measurements  discussed  previously  in  this  section. 

Applications  of  the  above  technique  to  a  propane-fueled  flame  indicated  a 
temperature  of  2241  +  30  K.  The  results  of  CARS  measurements  made  with  the 
ruby-based  CARS  system  are  listed  in  Table  3.  The  CARS  data  in  this  case 
indicated  a  temperature  of  2143.5  +  140  K.  The  large  variation  and  lower 
temperature  of  the  CARS  measurement  indicated  that  the  position  of  Lhe  flame 
was  somewhat  off  center  and  toward  the  more  turbulent  edges.  Another  measure¬ 
ment  was  conducted  on  the  burner  with  the  Nd:YAG-based  CARS  system  which  gave 
a  temperature  of  2220  K  +  50  K  for  the  same  flame  conditions  and  position  as 
the  sodium-line  reversal.  These  results  were  in  good  agreement  with  the 
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Temperature  Fit  of  CARS  Spectrum  Obtained  from  an  Acetylene-Fueled  Flame 


Table  3 


TEMPERATURE  FITS  TO  THE  PROPANE-FLAME  DATA 


Spectra  No. 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


Temperature 

00 

2009 . 7 

1853.7 

2389.9 

2252.9 

2177.4 

2208.2 

2125.8 

2059.9 

2010.6 

2086.4 

2185.4 

2315.9 

2167.3 
2117.0 

2215.5 
2309.8 

1876.4 

2220.5 


Dispersion 

(*K) 

31.86 

30.21 

30.67 

26.70 

26.95 

48.01 

24.07 

27.56 

34.81 

26.43 
41.78 
25.28 
30.72 
24.03 
25.99 
37.83 

41.43 

20.43 


Avg.  Temperature  =  2143.5  K  +  140  K 


reversal  measurements.  Temperature  profiles  which  were  obtained  on  the  flame 
are  shown  in  Fig.  31.  The  propane-fueled  flame  showed  the  same  basic  trends 
as  the  acetylene-fueled  flame — primarily  that  the  temperature  is  essentially 
flat  over  a  large  region  and  drops  off  sharply  at  the  edges.  Figure  31  shows 
a  typical  spectrum  and  temperature  fit  of  the  spectrum  for  these  data. 

3.7  REFERENCE-SCHEME  STUDIES 

In  order  to  obtain  single-shot  number  densities,  either  an  absolute  or  a 
relative  measurement  of  the  integrated  intensity  of  the  Raman  transition  must 
be  made.  The  former  is  extremely  difficult — if  not  impossible — in  most 
cases.  The  latter  is  performed  through  use  of  a  ratioing  technique  in  which 
the  majority  of  the  experimental  parameters  are  canceled.  In  an  effort  to 
discover  the  optimum  ratioing  arrangement,  a  preliminary  study  was  conducted 
on  various  signal-detection  schemes.  The  basic  experiment  consisted  of  a  ruby 
CARS  setup  in  which  the  sample  and  reference  legs  were  in  series  with  each 
other.  In  order  to  insure  that  the  CARS  signal  from  one  leg  would  not 
interfer  with  the  other,  special  blocking  filters  were  employed  which  passed 
the  ruby  and  dye  signals  but  blocked  the  CARS  signal.  Room  air  was  chosen  as 
the  gas  to  be  probed  for  both  the  reference  and  sample  paths.  Various  detec¬ 
tors  were  utilized  during  the  experiment  in  order  to  determine  the  best 
possible  arrangement. 

Results  of  this  experiment  are  given  in  Table  4.  A  ratioed  variation  of  10  - 
18%  was  encountered  depending  upon  the  particular  detectors  and/or  signals 
ratioed.  This  would  contribute  a  5  -  9%  error  in  the  measured  concentration. 
While  this  error  was  considered  to  be  relatively  small,  it  was  decided  that 
further  studies  into  lowering  these  errors  would  be  necessary .  An  important 
observation  was  made  during  the  incorporation  of  the  Nd:YAG  laser  into  the 
CARS  system — the  CARS  signal  fluctuations  increased  dramatically.  This 
observation  sparked  a  major  investigation  in  order  to  understand  this  phenome- 
num.  Figure  32  depicts  some  of  the  various  arrangements  that  were  employed 
during  the  testing.  In  Fig.  32(a)  a  parallel-path  configuration  was  used  in 
which  a  small  percentage  of  the  BOXCARS  beams  was  split  to  form  two  separate 
paths  in  which  the  CARS  signals  were  formed  by  two  independent  lenses.  This 
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Table  4 


RATIOING  ARRANGEMENT  AND  STANDARD  DEVIATIONS  OF  RATIOS 
OBTAINED  WITH  THE  RUBY  CARS  SYSTEM 


Standard  Deviation 


Ratio  Technique  _ of  Ratio 

OMA/OMA  11% 
OMA/Photomultiplier  16% 
OMA/PIN-Diode*  14% 
OMA/ (PIN-Diode)2  13% 
OMA/ (PIN-Diode) 3  12% 
Photomultiplier /Photomultiplier  10% 
Photomultiplier/PIN-Diode  18% 
Photomultiplier/ (PIN-Diode) ^  17% 
Photomultiplier/ (PIN-Diode) ^  17% 


The  PIN-Diode  monitored  the  pump  beam  during  the  experiment. 


I 


t 
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scheme  typically  showed  intensity  fluctuations  of  40  -  60%  in  each  leg. 
Ratioing  improved  the  variations  only  slightly.  The  slight  improvement  in 
the  signal  upon  ratioing  indicated  that  while  the  signals  were  correlated, 
the  degree  of  correlation  was  small.  The  small  degree  of  correlation  is 
believed  to  be  caused  by  the  inequality  of  the  sample  and  reference  paths. 
This  was  expected  to  be  true  especially  in  a  BOXCARS  arrangement,  primarily 
due  to  the  smaller  overlap  of  the  beams.  This  overlap  resulted  in  a  smaller 
region  for  CARS  production  in  which  small  path  differences  could  overlap 
different  regions  of  each  beam. 


The  basic  difference  between  the  Nd:YAG  and  the  ruby  laser  systems  which  were 
employed  for  these  studies  was  the  transverse  and  longitudinal  modes  of  the 
lasers.  The  ruby  laser  was  operated  under  single  transverse  (TEM^q)  and 
longitudinal  modes  (linewidth  of  0.035  cm  ^),  while  the  Nd:YAG  had  a  line- 
width  of  0.8  cm  ^  (~  160  longitudinal  modes)  and  operated  in  the  familiar 
donut  transverse  mode  typical  of  these  unstable  resonators.  Figure  33 
depicts  the  relative  mode  beating  of  the  longitudinal  modes  of  the  Nd:YAG  and 
dye  laser  as  recorded  with  a  100-ps  rise  time  photodiode.  The  modulation  of 
the  temporal  shape  of  the  Nd : YAG  is  quite  extensive.  The  modulation  shown  in 
the  dye  laser  is  due  to  the  off-axis  pumping  geometry  employed  with  the 
Nd:YAG.  The  modulation  is  somewhat  smoothed  due  to  the  averaging  effect 
(finite  response  time)  of  the  dye. 


The  CARS  intensity  which  is  governed  by  the  Fourier  transform  of  such  wave- 

32 

forms  as  shown  in  Fig.  33  is  given  by 

I  ,  =  Cl. I.  I.  fda  F(a  -  u>.)  [  d6  f  dy  |  X  (B  -  y)  |  2F(B  -  to  )  f  (y  -  u>  )  (39) 

jk£  j  k  ^  J  J  J  J  k 


where  I...  is  the  integrated  intensity  of  a  CARS  mode  generated  by  the  inter- 
J  K  V, 

action  of  two  pump  modes  with  spectra  I.F(a  -  to.)  and  I,F(8  -  m  )  and  one 

J  J  k  k 

Stokes  mode  I  f  (y  -  to  ) . 


The  quadratic  dependence  on  the  pump-mode  structure  and  linear  dependence  on 
probe-mode  structure  could  lead  to  small  changes  in  mode  structure  resulting 
in  large  changes  in  intensity.  Figure  34  demonstrates  the  effect  of  temporal 
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shape  (and  thus  mode  structure)  on  the  CARS  intensity;  it  also  indicates  that 
while  the  overall  Intensities  of  the  pump  and  probe  lasers  remained  approxi¬ 
mately  the  same  for  three  shots,  the  CARS  signal  fluctuated  tremendously. 

This  was  due  primarily  to  the  interaction  of  the  pump  and  probe  mode  structure. 

A  comparison  of  BOXCARS  and  collinear  CARS  temporal  responses  was  conducted 
using  the  arrangement  depicted  in  Fig.  32(b).  In  this  arrangement  the  two 
CARS  signals  are  produced  by  the  same  lens  and  except  for  the  interaction 
zones  being  different,  the  two  paths  can  be  considered  to  be  identical.  But 
as  indicated  by  Fig.  35,  the  variation  in  the  temporal  responses  of  the 
BOXCARS  and  collinear  CARS  indicates  that  even  though  the  same  lens  was  used 
to  produce  both  signals,  the  difference  in  interaction  zones  was  enough  to 
cause  a  large  variation  in  temporal  response  and  signal  levels.  This  further 
indicates  that  the  sample  and  reference  paths  should  be  constructed  in  a 
manner  to  obtain  minimum  path  difference. 

With  this  in  mind,  the  setup  in  Fig.  32(c)  was  devised  as  a  viable  arrange¬ 
ment.  The  BOXCARS  signals  produced  by  this  arrangement  were  formed  by  the 
same  lens,  dispersed  by  the  same  monochromator,  and  detected  by  matching 
photomultipliers.  Typical  results  achieved  using  this  arrangement  indicated 
variations  in  both  legs  on  the  order  of  50  -  70%,  but  with  a  reduction  in  the 
overall  modulation  by  ratioing  of  at  least  a  factor  of  3.  While  this  is  a 
significant  reduction  in  the  variation,  it  still  represented  a  large  uncer¬ 
tainty  in  a  concentration  measurement.  The  paths  still  were  not  exact  enough 
for  operation  of  a  laser  with  multi-model  structure.  One  possibility  for 
further  reducing  this  variation  was  operating  the  Nd:YAG  in  a  single  mode, 
which  had  been  used  with  good  results  in  the  ruby  system. 

An  Electronic  Line-Narrowing  (ELN)  device  was  used  for  operating  the  Nd:YAG 

in  single  mode.  This  device  acts  to  force  the  electro-optical  Q-switch  to 

behave  as  a  lossy  Q-switch,  which  allows  build-up  of  the  laser  pulse  to  occur 

over  a  long  period  of  time.  This  slow  build-up  allows  the  intracavity  etalon 

12 

to  be  more  efficient  in  discriminating  against  unwanted  modes.  With  the 
ELN,  deviations  of  the  CARS  signal  of  ~  7%  were  obtained  when  ratioing 
resonant  and  utilizing  the  configuration  of  Fig.  32(c).  This  experimental 
technique  reduced  to  less  than  ~  5%  the  uncertainties  of  concentration  measure¬ 
ments  due  to  laser-power  and  mode-structure  fluctuations. 
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3.8  SIMULTANEOUS  MEASUREMENTS 


In  the  past,  concentration  measurements  have  been  made  by  first  determining 

the  temperature  at  a  given  location  either  by  scanning  or  by  a  single  shot 

and  then  determining  the  concentration  by  integration  using  a  photomulti- 
33 

plier.  The  net  effect  has  been  that  temperature  and  concentration  determi¬ 
nations  were  made  at  different  times.  The  magnitude  of  the  error  in  the 
concentration  measurement  due  to  the  variation  of  the  temperature  can  be 
assessed  by  the  following:  The  mole  fraction  on  nitrogen  (xN  ),  as  deter- 
mined  by  a  CARS  measurement,  is  given  by 


=  0.7808 


where  T  =  T 

2  flame 

I^  =  intensity  at  T2 

T^  =  calibration  temperature 

1^  =  intensity  at  T^ 

F  =  temperature-dependent  correlation  factor 
N2 


The  term  Fj^  is  a  temperature-dependent  quantum-mechanical  function  which 

considers  the  ideal  gas  law,  number  densities,  electron-spin  degeneracies, 

„  „  34 

and  number  of  rotational  levels.  F^  ranges  from  1.0  at  =  300  K  to  2.0 
at  T2  =  2100  K.  Thus,  for  small  variances  in  T^ 


r 


r 


From  Eq.  (40)  the  probable  error  in  xj^  from  probable  errors  in  T^,  T^,  I^, 
and  1^  (assuming  no  intensity  losses  due  to  turbulence  effects)  is  given  by 
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(42) 


Combining  Eqs.  (40)  and  (42),  performing  the  differentiation,  and  simplifying 
yields 


2  2  2  2 

E  +  E  +  1/4  E  +  1/4  E 

1  2  L2  1 


1/2 


(43) 


where  E  is  the  relative  error  of  the  specific  parameter. 


From  this  relationship,  it  can  be  seen  that  the  major  contributors  to 
are  ET^  and  Thus,  the  associated  error  in  the  temperature  determination 

due  to  instrumental  factors  as  well  as  measurement  uncertainties  has  a 
tremendous  effect  upon  the  precision  of  the  concentration  measurement. 
Simultaneous  determination  of  temperature  and  number  density  assures  that  the 
relative  error  introduced  due  to  turbulence  during  time-averaged  measurements 
will  not  enter  into  the  concentration  measurement.  Simultaneous  measurements 
are  possible  because  the  general  shape  of  the  Q-branch  profile  can  be  used  to 
determine  the  temperature,  while  the  integrated  area  of  the  Q-branch  can  be 
used  to  determine  the  concentration. 


In  order  to  make  simultaneous  measurements,  the  sample  and  reference  signals 
which  are  produced  must  be  simultaneously  recorded  by  either  a  photomultiplier- 
OMA  combination  or  by  a  TN1710-DARSS  detector.  The  DARSS  detector  was  pri¬ 
marily  used  for  these  measurements  due  to  the  simple  detection  and  electronics 
scheme  it  offers. 


! 
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The  major  problem  encountered  with  the  DARSS  unit  was  the  single  dimensionality 
of  the  reticon  detector.  Previously  the  simultaneous  recording  of  a  sample 
and  reference  signal  required  that  each  signal  be  placed  on  a  separate  half  of 
a  two-dimensional  detector  such  as  the  PAR  1205  OMA.  However,  in  the  case  of 
the  DARSS,  a  different  arrangement  (shown  in  Fig.  36)  had  to  be  employed.  The 
two  beams  to  be  input  into  the  spectrometer  were  crossed  slightly  more  than 
one  focal  length  in  front  of  the  collection  lens.  This  allowed  the  two 
signals  to  run  parallel  and  also  focus  at  the  required  position  on  the  entrance 
slits.  The  slits  were  then  rotated  90°  from  their  customary  vertical  position 
into  a  horizontal  position  so  that  the  signals  could  be  placed  side  by  side 
and  simultaneously  recorded  by  the  linear  reticon  detector.  A  typical  display 
taken  during  the  simultaneous  measurements  is  shown  in  Fig.  37.  Data  were 
gathered  and  stored  in  the  Spectral  Data  Memory  (SDM)  of  the  DARSS  unit.  Once 
filled,  the  SDM  was  dumped  onto  a  floppy  disk  and  accumulation  continued.  All 
data  were  taken  under  ELN  operation.  The  timing  diagram  for  this  mode  of 
operation  is  discussed  in  Section  3.1. 

Simultaneous  measurements  of  temperature  and  concentration  were  obtained  on  a 
premixed  propane  fueled  flame,  and  the  results  are  shown  in  Figs.  38  and  39. 

in  the  air  was  used  as  the  reference  medium  in  these  concentration  studies. 
Figure  38  is  a  horizontal  profile  of  the  flame  starting  from  the  center 
position  (0  mm)  and  moving  outward.  The  figure  indicates  that  toward  the  edge 
of  the  flame  the  temperature  decreased  while  the  ^  concentration  increased. 

The  concentration  at  the  maximum  temperature  2255  K  of  both  profiles  was  found 
to  be  -  70%,  which  is  comparable  to  a  calculated  temperature  under  stoichio¬ 
metric  conditions  for  a  propane-air  premixed  flame  of  2275  K  and  a  concentra¬ 
tion  of  71. 5%  . 

Two  extremely  important  parameters  used  in  combustion  research  are  fuel-to- 

35 

air  ratio  and  combustion  efficiency.  Fuel-to-air  ratios  are  typically 
determined  by  gas-sampling  measurements  of  CO^,  CO,  NO^,  and  the  total 

hydrocarbon  concentration  (THC),  subject  to  the  procedures  set  forth  in 

36  * 

Aerospace  Recommended  Practice  1256.  Combustion  efficiencies  can  be 
determined  in  a  number  of  ways.  In  ramjet  combustor  development,  combustion 
efficiencies  are  calculated  from  thrust  or  temperature  measurements.  For 
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Figure  36.  Optical  Arrangement  for  Simultaneous  Recording 
of  a  Sample  and  Reference  Signal. 


Figure  32  Simultaneous  Measurement  of  Sample  and  Reference 
Signals  by  the  Tracor-Northern  DARSS. 
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Figure  38.  Horizontal  Profiles  of  the  Propane-Fueled  Flame 
from  the  Centerline  (0  mm)  of  the  Flame. 
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Figure  39.  Vertical  Profiles  of  the  Propane-Fueled  Flame. 
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gas-turbine  development,  combustion  efficiencies  normally  are  determined 

37 

from  emission  measurements  of  CO,  and  THC. 

The  CARS  technique  is  capable  of  making  ^  and  0^  concentration  measurements 

along  with  temperature  determinations  in  combusting  environments.  Assuming 

that  the  THC  levels  are  very  low,  the  following  equations  have  been  derived 

36 

for  the  determination  of  fuel-to-air  ratios  and  combustion  efficiencies: 


FA 


0.79  Mf 
(1  +  1/2  a)M 

a 


where 


a  =  hydrogen-atom-to-carbon-atom  ratios  in  the  fuel 

=  molecular  weight  of  fuel  =  12  +  a 

Ma  =  molecular  weight  of  air  =  28.9 

FA  =  fuel-to-air  ratio 

Xn  =  mole  fraction  of  0„ 

U2 

X„  =  mole  fraction  of  N0 
N_  t 


(«4) 


The  combustion  efficiency  is  given  by 


n  = 


T  -  T 

flame  inlet 


AT 


comb ,  ideal 


(45) 


where  AT  ,  , ,  ,  =  ideal  gas  temperature  calculated  from  the  fuel- 

comb  ,  ideal 


flame 


to-air  ratio 

=  measured  local  flame  temperature 


Tinlet  =  temperature  of  the  inlet  gases. 
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As  indicated,  the  above  equations  dictate  that  N^  and  0„  measurements  and 

2. 

flame-temperature  determinations  be  made  simultaneously  at  a  given  location. 
However,  in  the  case  of  a  premixed  flame  in  which  the  concentration  of  0^  is 
extremely  low  due  to  combustion,  only  the  N  concentration  and  temperature 
must  be  known  in  order  to  determine  the  fuel-to-air  ratios.  A  study  was 
conducted  using  the  premixed  burner  whose  profiles  were  reported  above.  The 
propane  fuel  flow  was  varied  and  the  simultaneous  temperatures  and  concentra¬ 
tions  determined  by  the  simultaneous  CARS  technique.  The  results  of  this 
study  are  shown  in  Fig.  40.  A*  indicated,  the  temperature  initially  is  low, 
maximizes,  and  then  gradually  decreases.  The  nitrogen  concentration,  however, 
shows  a  steady  decline  as  the  fuel  is  increased.  Given  the  temperature  and 
concentration,  the  fuel-to-air  ratios  can  be  determined  from  Eq.  (44).  A 
direct  comparison  of  the  calculated  fuel-to-air  ratios  could  not  be  made  with 
the  existing  flow  meters  since  they  were  calibrated  only  for  acetylene  flow 

and  not  propane.  However,  an  indirect  comparison  can  be  made  with  calculated 

38 

results  from  the  modified  version  of  the  NASA  ODE  program  for  calculating 
ideal  flame  temperatures,  fuel-to-air  ratios,  equivalence  ratios,  average 
molecular  weights  of  the  combustion  products,  and  equilibrium  concentrations 
of  CO^,  CO,  H^O,  H^ ,  0^  and  for  constant  pressure  combustion  for  a  specific 
fuel.  The  experimental  data  presented  in  Fig.  40  were  used  to  calculate  the 
fuel-to-air  ratio,  which  then  was  used  in  the  ODE  program  to  calculate  the 
ideal  flame  temperature  and  nitrogen  concentration.  The  results  of  this 
comparison  are  listed  in  Table  5.  The  results  agree  quite  well  with  the 
calculated  results,  which  indicates  that  this  is  a  viable  method  of  deter¬ 
mining  the  fuel-to-air  ratios. 

3.9  BACKGROUND- SUPPRESS ION  MEASUREMENTS  IN  PREMIXED  PROPANE  FLAME 


The  major  limitation  to  the  sensitivity  of  the  CARS  technique  is  the  inter¬ 
ference  effects  between  the  resonant  and  nonresonant  contributions  to  the 

third-order  susceptibility.  Several  methods  have  been  proposed  for  reduction 
39-42 

of  this  background.  Most  are  experimentally  complex  because  of  the  need 

for  three  lasers  or  the  need  for  excellent  polarization  characteristics; 
and,  in  some  cases,  only  one-order-of-magnitude  rejection  of  the  background  is 
produced.  Of  the  present  techniques  proposed,  polarization  cancellation 
is  considered  to  be  the  best  method  for  background  suppression  because  of  its 
relative  ease  of  application  and  complete  rejection  of  the  background.  The 
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Table  5 


COMPARISON  OF  TEMPERATURES  AND  NUMBER  DENSITIES  OBTAINED 
FROM  CARS  MEASUREMENTS  WITH  THOSE  PREDICTED 
BY  IDEAL  FLAME  CALCULATIONS  FROM  THE 
MEASURED  FUEL-TO-AIR  RATIOS 
(M  -  MEASURED,  C  -  CALCULATED) 


T 

m 

(%N,) 
l  m 

(F/A) 

m 

X 

c 

(%N0) 

2065 

74.0  +  4.3 

0.0595 

2209 

72.8 

2300 

71.0  +  6.5 

0.0701 

2265.7 

70.9 

2122 

70.5  +  7 

0.0718 

2250.7 

70.4 

2000 

66.7  +  10 

0.0858 

2087.4 

66.7 

2000 

68.3  +  5 

0.0796 

2165.0 

68.4 

1900 

64.4  +  4 

0.095 

1981.2 

64.3 

1850 

62.7  +  9 

0.1022 

1901.0 

62.8 
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basic  polarization  technique  for  a  two-laser  collinear  CARS  system  was 
reported  by  Rahn,  tit  al^.  In  this  setup  the  polarization  of  the  analyzer 
and  that  of  the  probe  beams  was  ~  60°  with  respect  to  the  pump,  as  shown  in 
Fig.  41(a).  This  arrangement  allowed  for  complete  cancellation  of  the  back¬ 
ground  with  a  maximum  of  resonant  signal  production. 

In  BOXCARS,  some  flexibility  is  afforded  by  the  availability  of  two  spatially 
distinct  laser  beams  at  which  can  have  different  polarizations.  Figure 
41(b)  shows  the  polarization  scheme  for  BOXCARS.  In  this  instance  the 
analyzer  angle  0  is  given  by 


0  =  tan 


3  cos8cos4>  +  siniftsinQ 
cos0sin<f>  +  sin0cos(j> 


where  the  CARS  intensity  is  proportional  to 

2 

I  oc  |  4  cos0cos<(>cos0  +  2  sin0sin<j>cos0  -  2  sin0cos<|>sin0 1  (47) 

Computer  solutions  of  Eqs.  (46)  and  (47)  as  a  function  of  0  for  the  maximum 
intensity  are  given  in  Table  6  and  plotted  in  Fig.  42.  As  indicated,  a 
minimum  loss  of  16  is  observed  with  background  suppression  as  compared  to  that 
of  the  parallel-polarization  arrangement  for  all  beams  which  is  normally  used. 

It  was  assumed  that  since  the  fuel  mixture  employed  by  the  burner  was  premixed, 

the  concentration  of  oxygen  inside  the  flame  was  very  low.  In  an  effort  to 

examine  oxygen,  rhodamine  610  was  incorporated  and  the  concentration  tuned  to 

the  correct  region.  Inside  the  propane-fueled  flame,  only  a  nonresonant 

background  could  be  seen;  but  as  the  burner  was  translated  to  the  edge  of  the 

flame,  the  interference  of  resonant  oxygen  and  nonresonant  background  was 

encountered.  Figure  43(a)  shows  the  spectrum  obtained  on  the  outer  edge  of 

the  flame.  In  order  to  clean  up  the  spectrum  and  to  permit  temperature 

determination,  background  cancellation  employing  the  technique  of  Rahn, 

43 

et  al.,  was  utilized.  The  results  are  shown  in  Fig.  43(b).  A  temperature 
of  1500  K  was  fit  to  this  spectrum.  While  no  attempt  was  made  to  obtain  the 


37 


Table  6 


COMPUTER  DETERMINATIONS  OF  BOXCARS  INTENSITY 
AND  ANGLES  OF  POLARIZATION  FOR 
BACKGROUND  SUPPRESSION 


_e_ 

f 

8 

I 

0 

0.0 

0.0 

16.00* 

0 

60.0 

60.0 

1.00 

5 

62.5 

57.5 

1.00 

10 

65.0 

55.0 

1.00 

15 

70.0 

51.0 

0.995 

20 

70.0 

52.0 

0.999 

25 

75.0 

48.5 

0.999 

30 

77.5 

47.5 

0.999 

35 

80.0 

47.6 

0.998 

40 

85.0 

45.7 

0.999 

45 

90.0 

45.0 

1.00 

50 

95.0 

46.05 

0.999 

55 

105.0 

45.32 

0.993 

60 

115.0 

60.0 

0.992 

65 

105.0 

72.39 

0.713 

70 

100.0 

76.91 

0.467 

75 

100.0 

83.9 

0.267 

80 

95.0 

84.67 

0.124 

85 

90.0 

85.0 

0.03 

90 

_ 

— 

0.00 

Parallel  polarization  employed  in  normal  CARS. 
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concentration  profile,  future  emphasis  should  be  placed  upon  mapping  the  flame 
in  the  lower-concentration  regions  in  order  to  acquire  the  lower  limit  of 
single-shot  detectivity. 

3.10  TURBULENCE  EFFECTS  ON  CARS  INTENSITIES  AND  A  SCHEME  FOR  THEIR  COMPENSATION 

The  CARS  intensity  depends  upon  tight  focusing  of  the  pump  and  Stokes  beams  at 
the  same  point  in  the  medium  under  investigation.  The  effect  of  turbulence 
upon  the  propagation  of  these  beams  is  to  cause  the  focal  volume  to  increase 

4  A 

with  a  corresponding  decrease  in  the  CARS  signal.  Density  fluctuations  of 
the  measured  species  concentration  in  a  turbulent  medium  also  cause  variation 
in  the  CARS  intensity  due  to  the  dependence  of  the  CARS  signal  upon  the  square 
of  the  number  density.  The  latter  variation,  which  is  the  desired  parameter 
of  the  measurement,  may  be  adversely  affected  by  the  effects  of  turbulence 
upon  the  beams.  Since  combustion  media  are  highly  turbulent,  it  is  important 
that  the  effects  of  turbulence  upon  the  operation  of  the  CARS  signal  be  crit¬ 
ically  evaluated  if  accurate  number  densities  are  to  be  measured. 

The  generation  of  CARS  occurs  in  a  gaseous  medium  when  two  intense  col linear 

•  electromagnetic  waves  interact  with  a  molecular  species.  The  higher- frequency 

wave,  E  (w  ) ,  which  is  called  the  pump  wave  and  the  lower- frequency  wave, 

E  (u>  -  <5),  generate  an  anti-Stokes  spectrum  at  2w  -  u>  +6  when  E  is  tuned 
s  s  p  s  s 

through  the  Stokes  spectrum  of  the  molecular  species  under  investigation.  The 
anti-Stokes  spectrum  can  also  be  generated  by  employing  a  broadband  source  for 

E  in  conjunction  with  the  narrow-band  pump  wave. ^  The  appropriate  term  in 

S  45 

the  third-order  polarization  which  is  responsible  for  this  effect  is 

P  (3)  (to  +  6)  =  3X(3)  («  ,  a)  »  +  6)  E^E  (ui  -  6)  (48) 

a  a  pps  pss 

(3) 

where  x  is  the  corresponding  third-order  susceptibility,  is  the  center 

frequency  of  the  Stokes  line,  6  is  the  deviation  from  w  ,  and  w  is  the  center 

45  s  a 

frequency  of  the  anti-Stokes  line  given  by  2w  -  wg.  Similar  third- order 
polarizations  can  be  written  at  the  pump  and  Stokes  frequencies.  However,  for 


9  ? 


T 


weak  interactions  the  amplitude  of  the  pump  and  Stokes  waves  are  not  signifi¬ 
cantly  affected  by  the  medium.  Thus,  only  the  electromagnetic -wave  equation 
at  the  anti-Stokes  frequency  must  be  considered  in  the  analysis.  This  wave 
equation  is 


_1_ 

2 

c 


+  4  IT  (P 

a 


(1) 


+  P 


(3). 


(49) 


where  E  is  the  anti-Stokes  field  and  P  ^  is  the  linear  polarization.  For 
a  a 

plane  waves,  the  solution  of  the  above  equation  yields  for  the  anti-Stokes 
intensity 


I  (co  +  6,  Z) 
a  a 


3X 


(3)  |  2 


2 


I 

P 


I  (u> 
s  s 


2 

sin  (ZAk/2) 
(Ak/2)2 


(50) 


where 


Ak 


2u  n  to  n  w  n 
p  p  s  s  a  a 

c  c  c 


and  n  and  n  are  the  indices  of  refraction.  For  most  applications  ZAk/2  is 
a  s 

small;  thus,  Eq.  (50)  reduces  to 


T  (w  +  fi,  Z) 
a  a 


|3x(3)|2Z2I^>s 


A) 


(51) 


To  obtain  spatial  resolution,  focused  beams  are  used  and  most  of  the  CARS 
intensity  is  generated  near  the  focal  volume.  For  this  case,  the  solution  for 
E  may  be  placed  in  the  integral  form^ 

a 


E  (r,  t)  =  1/c 
a 


« 


(3) 


(r\  t')  dr' 


(52) 


where  t'  *  t  -  |?  -  r'|/c.  As  can  be  seen  from  Eq.  (48),  the  evaluation  of 
this  integral  requires  an  explicit  expression  for  the  focused  pump  and  Stokes 

waves.  Although  approximate  analytical  expressions  for  E  and  E  are  avail- 

47  P  a 

able,  it  is  not  possible  to  carry  out  the  integration  to  yield  a  closed 
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form.  An  easier  approach  to  the  problem  utilizes  the  plane-wave  solution  for 

I  (cu  +  5)  at  the  focal  volume  and  simple  diffraction  theory  to  estimate  the 
a  a  48 

diameter  of  the  focal  spot.  This  approach  yields  the  approximate  solution 


P  (w  +5) 
a  a 


4iuo  a)  \  2 

-r) 


3X(3)|2  P2  P  (a) 

1  1  p  s  s 


6) 


(53) 


where  and  P^  are  the  corresponding  spectral  power  densities  and  P^  is  the 
power  of  the  monochromatic  pump. 


In  this  study,  number -density  measurements  are  made  by  using  a  Nd:YAG 
laser  and  a  broadband  Stokes  wave  from  a  broadband  dye  laser  which 
permits  the  simultaneous  generation  of  the  entire  Q-branch  of  the  molecular 
species  under  investigation.  The  duration  of  these  waves  is  approximately 
10  nsec,  and  a  photomultiplier  is  used  to  integrate  the  total  anti-Stokes 
power  in  the  Q-branch  for  each  pulse.  The  energy  in  a  pulse  may  be  expressed 


Int  " 


KN2F(T)  p2  p 
P  s 


eff 


(P.T) 


(54) 


if  the  contributions  due  to  nonresonant  species  are  negligible  and  Pg(ws  +  6) 
is  constant  over  the  Q-branch.  In  Eq.  (54),  N  is  the  number  density,  p  is  the 
total  pressure,  K  is  a  constant,  r  is  an  effective  linewidth  and  F(T)  is  a 
function  of  the  thermal  distribution  of  the  molecules  of  the  medium. 


A  complete  analysis  of  the  effect  of  turbulence  on  the  generation-  of  CARS 
requires  a  solution  of  Eq.  (52),  which  is  a  very  difficult  task.  However,  in  a 
turbulent  medium  the  focal  spot  increases  in  size  and  is  also  displaced  from 
the  nonturbulent  position  in  a  statistical  way  from  pulse  to  pulse.  The 
change  in  the  spot  size  should  have  the  greatest  effect  on  the  CARS  signal 
since  this  reduces  the  intensities  in  the  focal  volume.  For  this  case,  the 
plane -wave  solution  yields 
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at 


A  I 
t  a 


|3x<3,|2  22 


p  V 

i) 


P4 


PJ  (ws  -  6) 


(55) 


where  is  the  effective  area  of  the  focal  spot  and  Zg  is  the  corresponding 

distance  over  which  the  major  portion  of  the  CARS  signal  is  generated.  From 

geometrical  considerations,  it  is  reasonable  to  postulate  that  the  value  of  Zg 

increases  with  A  such  that 
e 

Z2  =  kA  (55) 

e  e 


where  k  is  a  constant.  Thus,  the  ratio  of  the  anti-Stokes  energy  for  a 
nonturbulent  medium  to  that  for  a  corresponding  turbulent  medium  can  be 
estimated  from  Eqs.  (55)  and  (56).  This  gives 

P  Pi.  A 

a  aj int  et 

P  ,  “  P  A  .  _  “  A 
at  atj  int  e 

where  the  subscript  t  is  used  for  parameters  corresponding  to  the  turbulent 
medium. 


et 


(57) 


Since  turbulence  affects  beam  propagation  by  inducing  random  fluctuations  in 
the  index  of  refraction  within  the  medium,  two  gases  having  a  large  difference 
in  their  refractive  indices  can  be  mixed  in  a  turbulent  jet  to  produce  a 
small  turbulent  field.  Thus,  to  investigate  the  effects  of  turbulence  upon 
beam  propagation  and  the  generation  of  CARS  signals,  the  turbulent  generator 
shown  in  Fig.  44  was  employed.  To  make  the  turbulence  field  produced  by  the 
generator  as  uniform  as  possible,  the  nozzles  of  the  generator  were  made  by 
photographically  producing  an  etching  pattern  on  copper-clad  epoxy  plates. 

The  copper  was  etched  away  in  a  pattern  which  formed  six  nozzles  on  each 
plate  when  the  plate  was  sandwiched  between  separating  plates.  The  nozzl 's 
were  fed  from  the  large  openings  in  the  plate  which,  in  turn,  were  fed  from 
the  header.  Alternating  the  orientations  of  the  plates  allowed  for  feeding 
even-  and  odd-numbered  nozzle  rows  with  two  different  gases  to  generate  a 
turbulent  mixing  medium. 


The  CARS  experimental  setup  for  studying  the  turbulence  effects  shown  in 
Fig.  45(a)  was  the  same  as  that  employed  for  the  simultaneous  temperature- 
species  concentration  measurements  discussed  in  Section  3.8,  with  the 
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Figure  45.  a)  Experimental  Setup  of  Turbulent  Studies 

b)  Setup  Employed  for  Broadband  and  Narrow-Band  Dye  Laser 


exception  that  a  collinear  rather  than  a  BOXCARS  arrangement  was  employed  for 
the  turbulence  studies.  The  turbulent  generator  in  these  studies  was  placed 
between  the  focusing  lens  and  its  focus  in  the  sample  leg  and  kept  well  away 
from  the  focal  region  where  the  majority  of  the  CARS  signal  is  produced. 

This  procedure  ensured  that  any  loss  in  the  CARS  signal  would  be  due  primarily 
to  turbulence  effects  and  not  to  number-density  changes  resulting  from  helium 
being  introduced  into  the  focal  volume. 

In  order  to  study  the  variation  in  spot  size  after  propagation  through  turbu¬ 
lence,  a  small  portion  of  the  pump  and  Stokes  beams  was  split  off  from  the 
sample  path  just  after  focus.  The  split  beams  were  then  coupled  onto  a  PAR 
1205  OMA  detector  which  allowed  recording  of  the  spread  of  the  focal  spot. 

Figure  46  shows  the  effects  of  turbulence  upon  the  spot  size  and  upon  the 
CARS  intensity  as  recorded  by  the  setup  described  above.  As  can  be  seen  from 
Fig.  46(a),  with  the  use  of  air  only  (and  no  helium  flow),  the  focal  spot  is 
sharp  and  well  defined;  the  resulting  CARS  signal  from  the  sample  path  is 
strong.  However,  when  the  helium  is  allowed  to  flow  (forming  a  refractive- 
index  gradient),  the  focal  spot  expands  drastically  with  a  corresponding 
decrease  in  the  CARS  intensity  in  the  sample  leg.  The  effect  of  introducing 
this  type  of  turbulence  in  the  beam  path,  in  this  case,  resulted  in  a  factor- 
of- eleven  reduction  in  intensity,  with  a  corresponding  330%  error  in  the 
concentration  measurement. 

Figure  47  shows  the  results  of  a  detailed  study  of  the  variation  of  the  CARS 
intensity  as  a  function  of  turbulent  generator  position.  Along  with  the 
relative  CARS  signal,  the  relative  beam  spread  as  measured  by  the  OMA  setup 
is  also  displayed.  The  x-axis  represents  the  distance  of  the  generator  from 
the  lens  for  a  16. 5-in. -focal-length  lens.  As  indicated  by  the  experimental 
curves,  the  relative  CARS  signal  decreases  as  the  distance  from  the  generator 
to  the  focal  region  of  the  beam  increases.  Inversely,  the  beam  spread  increases. 
A  quick  calculation  of  the  natural  log  of  the  relative  intensity  as  a  function 
of  the  natural  log  of  the  inverse  of  the  relative  beam  spread  indicates, 
however,  that  the  intensity  varies  not  as  the  square  of  the  inverse  of  the 
beam  spread  [see  Eq.  (59)]  but  simply  as  the  inverse  of  the  beam  spread. 

Figure  48  is  a  plot  of  the  experimentally  obtained  intensities  as  a  function 
of  the  intensities  predicted  by 
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and  as  Predicted  (x)  from  the  Measured  Beam 
“  1/D)  for  a  Horizontal  Scan  of  the  Turbulent 


(58) 


As  can  be  seen,  the  correlation  between  the  experimentally  observed  values 
and  those  predicted  by  Eq.  (58)  is  quite  good.  An  important  clue  concerning 
the  reason  for  this  result  was  obtained  when  it  was  noted  that  a  90-deg. 
rotation  of  the  OMA,  which  resulted  in  the  rectangular  target  being  perpen¬ 
dicular  to  the  beam  propagation,  gave  little  indication  that  the  beam  was 
spreading.  Imaging  the  beam  upon  the  laboratory  wall  allowed  the  visual 
assessment  that  the  beam  was  spreading  in  only  one  direction.  The  spread,  in 
fact,  was  perpendicular  to  the  flow  of  the  helium-air  turbulence,  indicating 
that  the  turbulence  was  not  completely  random  in  nature.  This  one-dimensional 
spreading  by  the  turbulent  field  results  in  a  linear  inverse  dependence  of 
the  CARS  intensity  upon  the  beam  spread,  as  indicated. 

Because  of  the  large  errors  that  can  be  introduced  into  the  concentration 
measurement,  it  is  imperative  that  some  means  of  compensating  for  this  turbj- 
lence  be  developed.  Several  such  schemes  have  been  proposed  and  were  con¬ 
sidered  for  this  study.  One  such  scheme  involves  taking  a  reference  after 
the  turbulence  has  been  introduced  in  the  beam,  thus  compensating  for  its 
effect.  However,  a  problem  with  this  method  is  that  while  the  CARS  signal  is 
affected  by  the  turbulence  through  which  the  beam  propagates  until  it  reaches 
the  focus,  it  is  unaffected  by  the  turbulence  which  occurs  after  the  focus. 
This,  however,  is  not  the  case  for  the  proposed  reference  because  turbulence 
which  occurs  after  the  focus  will  affect  its  CARS  intensity.  Such  a  reference 
would  tend  to  overcompensate  for  the  turbulence  since  the  reference  beams 
would  pass  through  twice  the  length  of  turbulence  that  the  sample  beams 
traversed  before  the  respective  CARS  signals  were  produced.  The  reference 
scheme  which  is  placed  before  or  parallel  to  the  sample  path  would  not  com¬ 
pensate  for  the  turbulence.  Therefore,  it  would  also  be,  at  best,  a  poor 
reference  for  a  turbulent  medium. 

What  was  considered  to  be  the  most  promising  approach  was  a  novel  reference 

49 

scheme  proposed  by  Ultee  and  reported  by  Shirley  for  background-free  CARS 
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spectra,  i.e.,  high  major  species  concentrations  or  low  concentrations  with 
the  background  cancelled.  This  scheme  involved  the  use  of  a  nonresonant 
signal  produced  at  the  measurement  location  as  an  In  situ  reference.  In  the 
above-cited  report,  it  was  proposed  that  a  two-component  dye  laser  be  employed 
to  produce  both  a  resonant  and  a  nonresonant  signal.  The  temperature  could 
be  determined  from  the  resonant,  and  the  nonresonant  could  serve  as  a  calibra¬ 
tion  source  for  the  concentration  measurement  since  it  is  a  measure  of  the 
total  number  of  molecules  present  in  the  measurement  volume.  It  also  was 
noted  that  such  a  measurement  approach  would  be  immune  to  medium  extinction 
and  refraction  effects  since  all  beams  would  be  similarly  affected,  thus 
making  it  an  ideal  method  for  turbulent  compensation.  The  article  fails  to 
mention  the  fact  that  a  power  reference  would  still  have  to  be  employed  due 
to  the  relative  power  fluctuations  which  are  inevitable  with  such  a  dye-laser 
system. 

Figure  45(b)  depicts  the  dye- laser  setup  used  for  producing  the  broadband  and 
narrow-band  Stokes  beams  needed  for  the  experiment.  A  grazing-angle  prism50 
was  used  for  beam  expansion,  thus  allowing  for  one-half  the  beam  to  be  inci¬ 
dent  upon  a  total  reflector  which  produces  the  broadband  Stokes  output  and 

O 

the  other  half  to  be  incident  upon  a  600-groove /mm  5000-A  blazed  grating 
which  accounted  for  the  narrow  nonresonant  Stokes  output.  The  output  of  the 
dye  laser  was  then  amplified  and  directed  to  the  turbulent  region,  as  shown 
in  Fig.  45(a). 

The  incorporation  of  the  additional  nonresonant  Stokes  output  results  in  four 
peaks  being  displayed  by  a  Tracor-Northern  1710-DARSS  multichannel  analyzer. ^ 
The  correlations  without  turbulent  flow  of  the  nonresonant  signals  in  the 
sample  and  reference  path  and  the  resonant  signals  in  the  sample  and  refer¬ 
ence  path  were  better  than  97%.  Figure  49  shows  the  results  obtained  with 
this  scheme  with  and  without  turbulent  flow.  Figure  49(b)  shows  the  display 
with  air  flow  only.  The  two  peaks  at  the  left  of  the  display  are  from  the 
reference  path  which  experienced  no  turbulence,  and  the  two  peaks  at  the 
right  are  from  the  sample  path  in  which  the  generator  was  located.  Figure 
49(a)  shows  the  results  when  allowing  helium  to  flow,  thus  forming  a  turbu¬ 
lent  mixing  region.  As  can  be  seen  the  peaks  in  the  sample  path  are  lower  in 
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Display  of  the  In  Situ  Reference  Scheme  a)  With  Flow  and  b)  Without  Flow 
N-Nonresonant  Peak,  R-Resonant  Peak. 


intensity  than  those  obtained  without  flow;  however,  both  the  nonresonant  and 
resonant  peaks  decreased  by  the  same  percentage.  By  normalizing  the  resonant 

t 

signal  to  the  nonresonant  in_  situ  reference,  one  can  compensate  for  the 
effect  of  turbulence.  Tables  7  and  8  show  the  experimental  results  that  were 
obtained  with  the  setup  described  above.  The  variation  of  the  relative 
signals  increased  from  ~  3%  to  -  9%  when  the  turbulence  was  activated,  but 
the  ratioed  results  were  ~  5%,  which  was  comparable  to  results  obtained 
without  turbulence.  Table  8  indicates  the  concentrations  that  would  have 
been  obtained  from  the  results  listed  in  Table  7,  with  and  without  compensa¬ 
tion  for  the  turbulence.  As  indicated  the  error  without  compensation  is 
quite  large,  while  the  compensated  concentration  is  in  good  agreement  with 
the  nonturbulent  result  for  ambient  air.  These  results  indicate  that  by 
using  this  method  for  turbulence  compensation,  accurate  intensity  information 
can  be  obtained  without  regard  to  the  severity  of  the  refractive-index  change 
or  beam  extinction.  Immediate  applications  for  this  type  of  method  lie  in 
concentration  measurements  in  turbulent  and  particle-laden  environments. 

This  method,  although  demonstrated  only  for  collinear  CARS,  should  also  be 
easily  applicable  to  BOXCAR-type  arrangements. 
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Table  7 


RESULTS  OF  TWO  SEPARATE  EXPERIMENTS  EMPLOYING  THE  IN  SITU 
TURBULENCE  COMPENSATION  TECHNIQUE 


No  Flow*+ 
Resonant 


Nonresonant 


Ratio  (R/N) 


Table  8 


APPARENT 


N2  CONCENTRATIONS  WITH  AND  WITHOUT  TURBULENCE  COMPENSATION 
~ (ASSUMING  78%  N2  CONCENTRATION  IN  AMBIENT  AIR) 


Without  Compensation 
(%) 

With  Compensation 
(%) 

60.82 

77.69 

68.94 

81.18 

65.05 

80.46 

66.85 

81.24 

67.76 

76.47 

68.00 

78.82 

x  =  66.24 

x  =  79.43 

0  =  2.96 

a  =  1.83 

+  4.47% 

+  2.30% 
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Section  4 


HARDENED-CARS-SYSTEM  DEVELOPMENT 

To  extend  the  study  of  the  CARS  diagnostic  techniques  developed  during  this 

4 

study  and  earlier  work  to  large-scale  practical  combustion  environments,  an 
environmentally  hardened  CARS  optical  system  was  required.  The  objective  of 
this  development  was  to  produce  a  CARS  system  which  could  obtain  temperature 
and  major-species-concentration  information  while  operating  in  the  relatively 
uncontrolled  environment  of  a  combustor  test  tunnel.  Following  is  a  summary 
of  the  system  design  considerations  and  development  of  the  CARS  diagnostic 
capabilities . 

In  designing  a  hardened  CARS  system,  several  environmental  problems  were 
anticipated.  Among  these  were  optical  problems  such  as  nonresonant  CARS 
generation  in  combustor  windows,  which  would  limit  the  detection  of  minority 
species,  and  mechanical  problems  connected  with  acoustic  and  seismic  vibrations 
and  thermal-gradient  stresses.  These  problems  were  attacked  through  the 
system  design.  This  design  included:  securely  bolting  all  components  to  a 
common  optical  table;  maintaining  the  optical  axis  as  close  to  the  table  as 
possible  to  reduce  bending  moments  in  optics  mounts;  merging  pump  and  Stokes 
beams  as  soon  as  possible  to  reduce  effects  of  relative  motion;  enclosing  the 
most  sensitive  components  and  the  system  as  a  whole  in  insulated  housings  to 
isolate  them  from  the  environment;  and  expanding  the  probing  beams  to  reduce 
CARS  generation  along  optical  paths  and  in  windows.  The  system  was  designed 
to  be  remotely  controlled  and  portable,  yet  rugged  enough  to  survive  in  the 
testing  environment  and  in  transport. 

Preparations  for  the  first  measurement  program  in  the  AFAPL  combustion 
tunnel  began  with  the  construction  of  the  initial  system  design  shown  in 
Fig.  50.  Temperature  measurements  required  spectral  scanning  of  both  sample 
and  reference  signals,  but  this  capability  was  not  achieved  in  this  system 
because  of  the  difficulty  involved  in  scanning  the  separate  sample  and 
reference  spectrometers  simultaneously.  Therefore,  the  initial  system,  which 
was  used  only  in  the  laboratory,  was  capable  of  performing  only  CARS  number- 
density  measurements. 
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Since  construction  of  the  first  system  went  quite  rapidly,  there  was  sufficient 
time  before  beginning  the  tunnel  experiments  to  modify  the  design  to  incorpo¬ 
rate  a  single  spectrometer  to  be  used  to  analyze  both  sample  and  reference 
signals.  This  design,  shown  in  Fig.  51,  uses  two  photomultipliers  for  detec¬ 
tion  of  the  sample  and  reference  signals  with  an  optical  multichannel  analyzer 
(OMA)  as  an  alternate  sample-detection  system.  The  two  photomultipliers 
mounted  on  a  common  spectrometer  provided  a  scanned  temperature  and  single¬ 
shot  concentration  capability,  while  the  OMA  was  used  to  facilitate  tuning  of 
the  Stokes  component  of  the  CARS  process.  The  primary  reason  for  the  incorpo¬ 
ration  of  the  OMA  was  to  provide  a  means  of  obtaining  a  complete  CARS  spectrum 
in  a  single  laser  pulse,  thus  allowing  the  possibility  of  single-shot  tempera¬ 
ture  determinations.  Once  this  design  was  constructed, the  system  was  tested 
in  the  laboratory  by  producing  CARS  signals  from  0^,  CO^,  CH^,  and  1^0 
molecules.  After  initial  testing  the  system  was  transported  to  the  combustor 
test  facility  to  begin  the  first  in  a  series  of  combustor  tests. 

The  optical  system  for  the  initial  test  series  is  shown  in  position  in  the 
combustion  facility  in  Fig.  52.  Also  shown  in  this  figure  is  the  optics- 
table  support  system  which  provides  three  degrees  of  movement  about  the 
combust ion -tunnel  axis.  The  CARS  electronics  control  console,  laser  cooler, 
and  power  supplies  are  shown  as  positioned  in  the  adjoining  control  room  in 
Fig.  53.  Preliminary  tests  in  the  three-month  program  were  designed  to 
exercise  the  system's  capabilities  in  varying  flame  conditions  to  determine 
its  ability  to  operate  in  the  combustor  environment.  The  results  indicated 
that  the  CARS  system  suffered  no  apparent  environmentally  induced  ill 
effects.  The  scope  of  the  test  program  was  then  extended  to  include  probing 
various  parts  of  the  flame  structure  in  an  attempt  to  identify  regions  where 
measurement  difficulties  such  as  background  luminescence  or  laser -induced 
radiation  from  soot  might  be  encountered.  This  study  was  performed  with  the 
aid  of  the  OMA  which  was  used  to  obtain  single-shot  and  0^  spectra  in  the 
flames  produced  by  several  fuel  types  and  flow  configurations.  A  slight 
difficulty  in  obtaining  "clean",  background-free  spectra  was  encountered  in 
the  extremely  luminous  portions  of  certain  flames.  This  problem  was  eliminated 
by  operating  the  OMA  in  a  1-ms  gated  mode.  In  all  cases  after  employing  the 
gated  OMA,  these  spectra  were  easily  obtainable,  indicating  complete  access  to 
combustion  regions  for  diagnostics  purposes. 
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Figure  51.  Second-Opt ical-System  Des 
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Scanned  temperatures  and  Nj  and  number-density  data  were  then  accumulated 
to  produce  profiles  of  these  parameters  along  the  combustor  centerline  for  a 
well-defined  flame  produced  by  a  fuel  nozzle.  Combustion-region  temperatures 
were  determined  using  the  spectrum  of  which  was  generated  by  stepping  the 
spectrometer  in  1-X  intervals  through  successive  anti-Stokes  spectra.  One 
such  spectrum  is  represented  by  the  +  signs  in  Fig.  54.  Each  of  these  points 
represents  the  average  of  30  consecutive  measurements  normalized  with  respect 
to  the  spectral-energy  density  of  the  reference  anti-Stokes  beam.  The 
resulting  composite  spectrum  was  then  compared  with  computer-generated 
spectra  at  various  temperatures  until  the  "best-fit"  temperature  was  found. 

The  solid  overlay  in  Fig.  54  shows  such  a  computer-generated  CARS  spectrum 
at  a  temperature  of  1700  K.  Since  this  spectral  scanning  and  overlaying 
technique  involves  various  time-averaging  effects,  the  temperatures  obtained 
are  necessarily  representative  of  a  "weighted-average"  value. 

Species  concentrations  were  determined  through  single-pulse  integrated  CARS 
intensity  measurements  involving  the  entire  Q-branch  of  the  ^  and  0^ 
constituents  of  the  combustion  gases.  The  data  from  these  measurements  were 
reduced  in  conjunction  with  the  "weighted-average"  temperatures  to  obtain 
percent-by-volume  concentrations.  The  temperature  and  concentration  profiles 
obtained  for  the  axial  test  positions  are  shown  in  Fig.  55.  Although  not 
indicated  on  this  plot,  the  minimum  detectable  concentration  capability  of 
the  CARS  system  in  these  experiments  was  determined  to  be  ~  2%  for  and 
0.5%  for  C^.  A  more  complete  description  of  these  and  other  results  of  this 
test  program  is  contained  in  Ref.  5.^ 

After  the  first  series  of  tunnel  tests,  the  optical  system  was  modified  to 
allow  a  horizontal  flame-profiling  capability  in  addition  to  the  vertical  and 
longitudinal  freedom  already  available.  This  capability  was  incorporated  so 
that  during  the  following  two-month  testing  period,  horizontal  temperature 
and  concentration  profiles  could  be  measured. 

The  purpose  of  this  series  of  experiments  was  to  obtain  profiles  in  a  given 
flame  by  means  of  the  CARS  and  conventional  probe  techniques  to  allow  quanti¬ 
tative  comparison  of  these  measurements.  Figure  5b  contains  the  CARS-developed 
profiles  obtained  through  the  flame  diameter  at  a  given  axial  position.  A 
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Pt  13%-RhFt  thermocouple  provided  by  NASA  Lewis  Research  Center  was  used  to 
profile  the  propane  flame  along  the  same  diameter.  Comparison  of  the  CARS 
and  thermocouple-derived  temperature  profiles  is  shown  in  Fig.  57.  The 
comparison  of  CARS-derived  0^  concentrations  and  those  obtained  with  a  gas¬ 
sampling  probe  is  shown  in  Fig.  58.  A  more  detailed  description  of  the 
experimental  conditions  and  data-reduction  considerations  and  a  discussion 
of  the  results  are  given  in  Ref.  33.  As  can  be  seen  from  these  comparative 
measurements,  the  agreement  between  the  two  methods  of  measurement  is  quite 
good  near  the  flame  centerline.  Significant  differences  appear,  however, 
for  positions  greater  than  approximately  4  cm  for  the  centerline  of  the  14- 
cm-diam.  flame.  In  the  temperature  comparisons  of  Fig.  57,  part  of  this 
discrepancy  was  attributed  to  the  effects  of  improper  "weighting"  of  the 
observed  spectra  introduced  by  the  spectral  scanning  technique  employed  for 
temperature  measurements.  To  eliminate  this  effect  and  also  to  provide  a 
means  of  obtaining  true-time  averaged  temperature  information,  a  single-shot 
"instantaneous"  temperature  capability  was  needed. 

Initial  attempts  at  single-shot  temperature  determination  within  the  combustor 
resulted  in  average  temperatures  that  were  consistently  lower  than  those 
measured  earlier.  The  lower  temperatures  were  determined  to  be  an  effect  of 
the  spectral  shape  of  the  broadband  dye  laser.  Figure  59  shows  successive 
shots  of  a  nonresonant  reference  signal  and  indicates  two  characteristics  of 
the  dye  laser  output  which  influence  single-shot  spectra.  The  shot-to-shot 
profile  changes  visible  in  the  figure  affected  measurement  repeatability, 
while  the  relative- intensity  distribution  about  the  resonant  frequency  had  a 
"weighting"  effect  on  the  generated  CARS  spectrum.  This  effect  of  nonuniform 
"pumping"  around  the  resonant  peak  caused  by  improper  tuning  and  stability 
of  the  frequency  of  the  dye  laser  was  believed  to  be  partly  responsible  for 
the  lowered  average  temperatures. 

In  order  to  correct  a  sample  spectrum  obtained  with  something  other  than  a 
"flat-topped"  dye-laser  profile,  a  frequency  normalization  of  the  sample  wave¬ 
form  was  employed.  This  normalization  capability  was  accomplished  by  using  a 
dual-channel  OMA  detector  which  allowed  two  spectra  to  be  recorded  simulta¬ 
neously.  Thus,  a  sample  waveform  and  a  nonresonant  reference  were  recorded 
and  a  channel-by-channel  intensity  normalization  of  the  sample-to-reference 
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Figure  57 .  Comparison  of  CARS  (Solid)  and  Thermocouple 
(Dashed)  Temperatures. 
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Figure  59.  Shot-to-Shot  Frequency  Fluctuation  in  Dye  Laser. 
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signal  was  performed.  The  photographs  in  Fig.  60  demonstrate  this  dual¬ 
spectrum  capability.  Figure  60(a)  is  a  room-temperature  resonant  N,,  CARS 
spectrum.  Figure  60(b)  is  the  nonresonant  CARS  generated  in  a  glass  substrate 
during  the  same  10-ns  laser  pulse.  This  spectrum  indicates  the  relative 
frequency  profile  of  the  Stokes  laser  output.  The  absolute  frequency  align¬ 
ment  between  the  resonant  and  reference  spectra  can  be  seen  in  Fig.  60(c). 

This  display  represents  the  resonant  signal  minus  the  nonresonant  and  indicates 
the  relative  ease  with  which  the  proper  tuning  of  the  Stokes  laser  wavelength 
can  be  accomplished.  Proper  tuning  ensures  that  sufficient  Stokes  power 
exists  for  pumping  of  the  higher  rotational  and  vibrational  energy  levels  of 
the  resonant  molecule.  To  accomplish  the  spectrum  normalization  as  well  as  to 
speed  up  the  recording  and  storage  of  the  large  amounts  of  OMA  information, 
the  OMA  was  interfaced  to  the  9030A  calculator.  Several  basic-language 
computer  programs  were  created  to  collect,  record,  and  analyze  the  spectral 
information.  Figure  61(b)  displays  the  frequency-normalized  1700  K  spectrum 
obtained  from  the  resonant  and  nonresonant  signals  of  Fig.  61(a)  and  treated 
in  the  manner  described  above. 

At  this  point  the  third  and  final  system  design  changes  were  made.  Since  a 
single-shot  temperature  capability  was  available,  only  slight  system  modifi¬ 
cations  were  necessary  to  acquire  the  sample  and  reference  spectra  and  their 
integrated  intensities  from  the  two  photomultipliers  on  a  shot-to-shot  basis, 
thus  allowing  simultaneous  single-shot  temperature  and  concentration  informa¬ 
tion  to  be  obtained.  An  additional  modification — placing  the  reference  genera¬ 
tion  cell  after  the  sample  volume — was  made  to  provide  better  compensation  for 
and/or  demonstration  of  turbulence  and  Schlieren  effects  produced  in  the 
sample  beam  path.  To  date  no  ill  effects  from  these  disturbances  have  been 

observed.  The  resulting  hardened  CARS  optical  system  design  is  shown  in  Fig. 

52 

62.  An  equipment  manual  describing  the  system  configuration  and  operation 
has  been  prepared. 

Characterization  of  the  final  system  configuration  was  carried  out  in  the 
laboratory  with  special  emphasis  on  determining  the  quality  of  the  measured 
parameters.  Several  sets  of  data  were  accumulated  to  give  an  indication  of 
the  accuracy  and  precision  with  which  the  hardened  system  could  perform 
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Figure  61,  (a)  Simultaneous  Resonant  and  Nonresonant  N.?  CARS  Signals 

(b)  Frequency-Normalized  N2  Spectrum  Taken  at  1700  K. 
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and  C>2  concentration  determinations.  These  results,  wllich  are  summarized  in 
Tables  9  and  10,  were  obtained  by  measurements  made  at  room  temperature  and 
1  atm.  of  pressure  in  a  25-cm-long  *  2.5-cm-diam.  cell  containing  ^  and  O2 
gases  mixed  in  the  stated  proportions.  As  can  be  seen  from  this  information, 
the  accuracy  of  the  measurements  degrades  at  the  lowest  concentrations  of 
both  species,  because  of  nonresonant  background  contributions. 

To  characterize  the  simultaneous  single-shot  temperature  and  concentration 

capability  of  the  hardened  system,  a  high-temperature  sample  environment  was 

created  by  placing  a  1. 3-cm-diam.  x  15-cm-long  stainless-steel  tube  above  a 

propane-fueled  laboratory  burner.  A  Pt  10%-RhPt  thermocouple  was  inserted 

parallel  to  the  tube  axis  to  a  point  =  2  mm  below  the  focus  of  the  sample 

beam.  With  this  arrangement,  "oven"  temperatures  of  up  to  1000  K  could  be 

attained.  No  thermocouple-radiation  corrections  were  applied  in  these 

measurements  since  the  required  corrections  of  these  temperatures  were  believed 

to  be  lower  than  other  experimental  uncertainties.  Since  the  ends  of  the 

oven  extended  well  beyond  the  flame  zone,  it  was  assumed  that  no  combustion 

products  were  introduced  into  the  tube.  Thus,  the  mole  fraction  of  ^ 

within  the  tube  should  remain  constant  at  78.1%  in  atmosphere,  while  the 

_2 

number  density  and,  hence,  CARS  intensity  displays  a  T  dependence.  The 
results  of  these  measurements  are  summarized  in  Table  11. 

Another  parameter  of  the  hardened  CARS  system  which  was  characterized  was  the 
volume  of  sample  gas  which  is  active  in  the  production  of  CARS  signals  in  the 
given  collinear  optical  configuration.  The  profile  shown  in  Fig.  63  describes 
the  relative  intensity  of  resonant  CARS  generation  as  a  function  of  distance 
from  the  focal  point  in  the  sample  volume.  This  curve  was  generated  by 
observing  the  integrated  resonant  C02  CARS  intensity  produced  when  a  rectan¬ 
gular  "sheet"  of  C02  gas,  1  mm  in  thickness,  was  translated  along  the  sample 
optical  path.  These  data  indicate  that  if  the  concentration  of  the  species 
being  measured  remains  fairly  uniform — between  +3.5  cm  and  -3.5  cm  from  the 
sample  focus— then  *  99%  of  the  generated  CARS  signal  will  be  produced  within 
a  path  length  of  12  mm.  Thus,  the  active  sample  volume  can  be  approximated 
by  a  cylinder  12  ram  in  length  and  a  focal  diameter  which  has  been  calculated 
to  be  ~  0.03  mm. 
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RELATIVE  CARS  INTENSITY 


Section  5 


CONCLUSIONS 

Numerous  fundamental  studies  have  been  conducted  with  the  laboratory  CARS 
system  to  determine  the  ability  of  the  technique  to  measure  temperatures  and 
number  densities  in  flame  environments.  Major  emphasis  has  been  placed  upon 
making  single-shot  thermometry  and  density  determinations.  The  results  of 
single-shot  thermometry  studies  on  a  propane-fueled  burner  have  been  compared 
with  those  obtained  by  means  of  the  sodium-line-reversal  technique  with  good 
agreement.  The  simultaneous  temperature-number  density  results  reported 
here  showed  good  agreement  with  adiabatic  flame  calculations.  The  effects 
of  turbulence  upon  CARS  signal  production  have  been  examined,  and  a  method 
for  compensation  has  been  demonstrated.  Major  emphasis  for  the  future 
should  be  in  the  areas  of  broadening  the  range  of  quantifiable  species  for 
the  CARS  technique,  increasing  the  data-acquisition  rates  of  the  CARS 
instrument  to  permit  time  histories  of  temperature-number  density  variations 
to  be  obtained,* and  combining  the  CARS  technique  with  other  techniques  such 
as  Laser  Doppler  Velicometry  (LDV)  to  permit  determination  of  velocities  and 
thus  turbulence  characteristics  of  combustion  systems. 

The  efforts  in  the  hardened-CARS-system  development  portion  of  this  program 
have  culminated  in  an  optical  diagnostic  system  having  a  demonstrated  ability 
to  survive  and  to  perform  temperature  and  species  concentration  measurements 
in  a  practical  combustion  environment.  The  capability  of  obtaining  these 
important  parameters  simultaneously  and  within  a  single  shot  has  also  been 
demonstrated.  These  results  represent  a  very  significant  advancement  for 
combustion-diagnostic  techniques;  however,  several  areas  for  further  develop¬ 
ment  can  be  recommended.  Among  these,  refinements  in  the  development  of 
CARS  theory  are  required  for  more  accurate  prediction  of  experimental 
observations.  For  measurements  with  finer  spatial  resolution,  a  CARS 
technique  using  small-angle  CARS  or  BOXCARS  should  be  considered  to  replace 
the  collinear  configuration  used  in  the  present  system.  Also,  additional 
study  is  required  to  determine  the  precision  and  accuracy  with  which  CARS 
measurements  can  be  made  in  a  combustion  environment. 
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LISTING  OF  THE  CARS  TEMPERATURE 
FITrTNG  PROGRAM 
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